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THEORETICAL  STUDIES  ON  FREE  ELECTRON  LASERS 


I .  SUMMARY 

/ 

J 

/ 

At  the  beginning  of  the  contract  year  (August  1, 
1981  to  August  1,  1982)  we  undertook  analytic  and 
computational  study^of  certain  theoretical  questions 
relating  to  Free  Electron  Lasers  (FEL's).  Ao  originally 
conceived,'  the  study  was  divided  into  five  tasks: 

(l^)  A  proper  Hamiltonian  formulation  of  the 
nonlinear  equations  for  gain-expanded  FEL's  will  be  made 
and  Liouville's  theorem  applied  to  these  systems. 

< 2.  The  analytic  theory  for  gain-expanded  FEL's 
will  be  extended  to  include  variable  property  gain 
expanders  as  proposed  for  storage  rings. 

3.  The  results  of  Tasks  1  and  2  will  be  applied, 
numerically  if  necessary,  to  the  determination  of  whether 
a  steady  state  gain-expanded  FEL  is  compatible  with 
storage  ring  operation. 

4. '  The  gain  of  a  variable  pitch  wiggler  (as  a 
function  of  signal  amplitude  and  frequency)  will  be 
studied  using  a  single  particle  simulation  code. 

5i  A  fully  self-consistent  particle  and  EM  wave 
( 1-D)  code,  with  provisions  for  treating  micro  bunches 
and  providing  frequency  discrimination  in  order  to 
consider  Raman  instability  effects,  will  be  developed  to 


,  treat  the  startup  and  pulse  formation  problem  for 
variable  pitch  wigglers.^N 

The  first  three  concerned  gain-expanded  FEL's  and 
the  last  two,  FEL's  with  variable  parameter  wigglers. 

Tasks  1  and  2  have  essentially  been  completed.  A 
Hamiltonian  formulation  of  the  electron  equation  of 
motion  for  the  "thin  lens"  gain-expanded  FEL  has  been 
derived  (Appendix  A) .  This  leads  to  two  fundamental 
limitations  on  gain  expansion  schemes:  (a)  The  Manley 
Rowe  relations,  and  (b)  generalized  gain-spread  theorems. 
Having  established  these  limits.  Task  3  remains  to  be 
completed.  The  implications  of  these  limits  with  respect 
to  the  overall  efficiency  of  gain-expanded  FELs  in 
storage  rings  appear  to  be  somewhat  pessimistic. 

Considerably  more  time  and  effort  were  devoted 
to  Task  4  and  Task  5  than  had  been  planned.  Consequently, 
it  was  possible  to  carry  these  tasks  much  further  than 
anticipated.  The  linear  eigenmode  analysis  leading  to 
a  calculation  of  the  linear  gain  as  a  function  of  the  FEL 
parameters  has  been  completed.  The  results  for  the  limit 
of  very  short  electron  pulses  are  discussed  in  Appendix  B. 
By  very  short  electron  pulses  is  meant  pulses  with  pulse 
length  £  much  shorter  than  the  slippage  distance 
kwLw/kg  >  r**£,  where  r-  2kwLwAYr/Yr  and  AVYr  is 
the  fractional  change  in  the  resonant  energy  Yr  down 
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the  wiggler.  A  discussion  of  the  results  for  long 
electron  pulses  is  under  preparation  and  will  be 
available  in  the  very  near  future.  A  one-dimensional 
self-consistent  particle  and  electromagnetic  wave  code 
to  simulate  the  operation  of  an  FEL  with  variable 
parameter  wigglers  and  with  frequency  filtering  of  the 
electromagnetic  pulse  has  been  written.  The  details  of 
this  simulation  code  are  described  in  Appendix  C.  This 
code  enabled  us  to  study  not  only  the  startup  and  pulse 
formation  problem,  but  also  the  stability  of  the 
saturated  state.  The  results  are  discussed  in 
Appendix  C. 

Our  principal  results  are: 

1.  We  have  derived  analytic  results  for  linear 
pass  to  pass  growth  of  the  pulse  at  low  amplitude.  For 
very  short  pulses,  these  are  given  in  Appendix  B.  These 
results  have  been  numerically  verified  in  Appendix  C. 

2.  We  find  that  for  high  saturation  systems 
optical  discrimination  is  required  in  order  to  prevent 
nonlinear  signal  breakup.  Qualitative  agreement  with  the 
theory  of  sideband  instablities  is  obtained  for  the  width 
of  the  required  frequency  filter. 

3.  The  further  requirement  that  the  linear 
phase  of  growth  saturate  at  a  high  enough  level  to  ensure 
trapping  seems  to  be  satisfied  except  for  systems  which 
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seek  to  enhance  linear  gain  by  use  of  a  long  constant 
parameter  section. 

4.  With  proper  choice  of  frequency  discrimination 
(and  for  an  ultra-short  electron  pulse)  we  have  demon¬ 
strated  very  high  extraction  systems  (30%  efficiency  at 
ly)  which  grow  from  noise  to  a  stable  steady  state.  We 
derive  and  numerically  verify  a  criterion 

(1  -  r)  Ay 

>  1.28  x  10  - -  Amps 

Yr 

for  the  minimum  current  required  for  a  system  with 
reflection  coefficient  r. 

5.  The  limitation  on  beam  thermal  energy  spread 
Yth  is  determined  more  by  the  need  for  good  linear  gain 
(Yth/r55  <  l)  than  by  that  for  effective  electron 
trapping,  (7th/ris  <  2)  • 

The  ultra-short  pulse  was  modeled  by  a  6-function 
in  the  simulation.  Due  to  computer  time  limitations, 
this  limit  is  easier  to  simulate  since  the  number  of 
particle  orbits  which  must  be  followed  is  far  fewer  for  a 
single  electron  beamlet  than  for  the  many  beamlets 
required  to  represent  a  long  continuous  electron  pulse. 
This  limit  is  also  analytically  tractable.  However,  the 
basic  physics  of  the  electron-photon  interaction  is 


expected  to  be  qualitatively  the  same  for  long  electron 
pulses.  Preliminary  simulations  of  long  electron  pulses 
have  indeed  exhibited  a  behavior  qualitatively  similar  to 
that  of  short  pulses.  The  corresponding  criterion  for 
the  minimum  current,  however,  may  be  slightly  more 
severe. 


Gain-Expanded  FEL's 

In  constant  parameter  FEL's,  electrons  initially 
lose  energy  to  the  radiation  field,  but  in  so  doing 
experience  a  gradual  shift  in  relative  phase  which 
eventually  reverses  the  direction  of  energy  flow  and 
impairs  gain  of  the  PEL.  To  overcome  this  limitation. 
Smith,  et  al.1  have  outlined  a  scheme  for  a  gain-expanded 
FEL  in  which  the  sensitivity  of  gain  to  variations  in  the 
energy  of  the  exciting  electrons  is  reduced.  This  scheme 
uses  a  periodic  magnet  with  transverse  gradients  such 
that  the  electrons  which  lose  energy  move  to  different 
transverse  positions  and  continue  to  remain  in  phase  with 
the  transverse  electric  field  of  the  electromagnetic 
radiation.  The  reduced  sensitivity  to  energy  enhances 
the  prospects  for  operation  of  a  gain-expanded  FEL  in 
conjunction  with  a  storage  ring,  since  many  passes  of  the 
circulating  electrons  through  the  FEL  would  occur  before 
gain  is  significantly  affected  by  energy  spread. 


However,  transverse  betatron  oscillations  of  the  electron 
trajectories  are  simultaneousely  excited  and  may  negate 
the  effectiveness  of  gain  expansion  unless  the  oscilla¬ 
tion  amplitudes  are  kept  small.  For  steady  state  opera¬ 
tion  in  a  storage  ring,  it  is  desirable  to  maximize  the 
gain  while  minimizing  the  level  of  betatron  oscillations. 
In  this  mode,  a  steady  state  is  conceivable  in  which  the 
growth  in  betatron  oscillations  in  the  FEL  is  damped  by 
incoherent  synchrotron  radiation  in  the  storage  ring. 

The  important  issue  to  be  resolved  is  whether  a  steady 
state  can  be  extablished  with  a  reasonable  ratio  of 
energy  gain  in  the  laser  to  energy  loss  in  synchrotron 

radiation  <A£laser/A£sync)  . 

2 

Madey,  et  al.  have  considered  the  feasibility  of 
cancelling  transverse  betatron  oscillations  by  a  careful 
design  of  the  FEL  magnet.  Their  results  on  the  gain 
characteristics  of  gain-expanded  FEL's  with  excitation 
cancellation  were  obtained  from  numerical  integration  of 
the  equations  of  motion.  In  our  view,  this  procedure  did 
not  provide  a  satisfactory  theoretical  understanding  of 
the  effect  of  excitation  cancellation  on  gain. 

Excitation  cancellation  seemed  to  involve  detailed 
cancellation  of  "big"  terms  in  the  equations,  and  it  was 
not  clear  what  residual  effects  were  responsible  for 
gain. 
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In  order  to  gain  some  insight  regarding  the  gain 
characteristics  of  gain-expanded  FEL's,  we  have 
investigated  the  "thin  Lens"  gain  expander  as  recently 
proposed  by  Madey,  et  al.^  using  a  Hamiltonian  formulation 
of  the  electron  equations  of  motion.  The  thin  lens 
Hamiltonian  was  transformed  by  a  series  of  canonical 
transformations  to  an  action-angle  Hamiltonian  of  the' 
form 


K  =  k«j  J  + 


lco 

ee  s 


(t  +  {1/6h-1)z) 


h  ( J ,  y )  e 

n,m 


i  (nk  z  +  m<j> ) 


w 


where  J  is  the  phase  area  for  betatron  oscillations  and  <J> 
its  conjugate  angle,  Y  is  the  electron  energy  and  x-z  its 
conjugate  "momentum".  z  is  the  independent  variable  and 
e  is  proportional  to  the  laser  field. 

The  Liouville  equation  for  the  electron 
distribution  function  may  now  be  written  down, 
solutions  obtained  by  expanding  in  powers  of  e,  and  the 
changes  <  6 J  >  and  <  6y  y  after  passage  through  the 
thin  lens  wiggler  determined.  The  details  of  the 
analysis  are  discussed  in  Appendix  A. 

In  the  usual  mode  of  FEL  operation  where  one 
resonance  is  of  importance  (which  is  practical,  necessary 
and  desirable  to  get  cancellation  of  excitation  at  a 
definite  z  value)  the  Manley-Rowe  relation  holds 
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<5J> 


-  m 


<Sy> 

0) 

s 

where  m  ?  0.  Thus,  gain  is  proportional  to  betatron 
excitation. 

The  m  *  0  resonance  does  not  produce  betatron 
excitations,  but  since  this  depends  on  energy 
derivatives,  presumably  large  energy  spread  results. 

A  generalization  of  the  Madey  theorem  relating 
gain  to  spread  (without  making  the  resonant  approxima¬ 
tion;  that  is,  keeping  all  m  values)  can  readily  be 
derived 


1  J_ 

2  9J 


<(SJ,!>  *  I  £  <SJ5V> 


1  J_ 

2  3y 


<<SyO  +  1  J3  <6J6Y> 


These  relations  impose  severe  restrictions  on 
gain  expander  operation 

laser 
sync 


e 


where  <$£  is  the  maximum  beam  spread  allowed  in  the 
storage  ring.  At  the  present  time,  the  existing 
proposals  for  operating  a  gain-expanded  FEL  in  a  storage 
ring  do  not  seem  capable  of  yielding  overall  efficiencies 
of  more  than  a  few  percent. 

Variable  Parameter  FEL'S 

An  alternative  scheme  for  overcoming  some  of  the 

limitations  of  constant  parameter  FEL's  is  the  variable 

4 

parameter  wiggler  discussed  by  Kroll,  et  al.  In  this 
scheme,  a  significant  fraction  of  electrons  in  the  beam 
is  trapped  in  the  ponderomotive  potential  well  produced, 
by  the  combined  wiggler  and  radiation  fields.  The 
resonant  energy  of  the  FEL  is  reduced  by  suitable 
variations  of  the  wiggler  parameters  and  this  in  turn 
results  in  the  extraction  of  energy  from  the  trapped 
electrons.  This  reduction  in  energy  can  be  much  larger 
than  the  energy  which  can  be  transferred  in  constant 
parameter  wigglers .  Thus,  the  use  of  variable  parameter 
wigglers  to  produce  high  power  radiation  is  a  promising 
prospect.  The  primary  objectives  of  our  investigation  of 
FEL's  with  variable  parameter  wigglers  are  to  define  the 
area  in  parameter  space  where 

1.  The  operation  of  the  FEL  can  be  initiated  by 
growth  of  the  electromagnetic  pulse  from  noise  levels. 
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2.  The  nonlinear  saturation  of  pulse  growth 
occurs  with  the  trapping  of  a  significant  fraction  of 


beam  electrons. 

3.  The  finite  amplitude  saturated  pulse 
propagates  without  deterioration  due  to  trapped  particle 
sideband  instabilities. 

4.  The  energy  extraction  from  the  electron  beam 
is  optimized,  consistent  with  reasonable  growth  rates, 
effective  trapping,  and  stable  propagation. 

To  achieve  these  objectives,  we  have  investigated 
in  detail  the  operation  of  an  FEL  configuration  in  which 
a  combination  of  constant  and  variable  parameter  wigglers 
is  positioned  between  mirrors.  An  electromagnetic  pulse 
generated  in  the  wiggler  is  reflected  repeatedly  by  the 
mirrors,  and  multiple  passes  of  the  pulse  through  the 
wiggler  are  obtained.  Relativistic  electron  beamlets 
are  injected  into  the  wiggler  at  periodic  intervals 
so  that  on  each  forward  pass  of  the  pulse  through  the 
wiggler,  there  is  overlap  of  the  pulse  and  beamlets. 

The  pulse  grows  on  each  forward  pass.  Eventually,  the 
amplitude  of  the  pulse  becomes  large  enough  to  trap  a 
fraction  of  the  beam  electrons.  At  this  stage,  the  pulse 
is  vulnerable  to  distortions  due  to  the  onset  of  trapped 
electron  sideband  instabilities.  In  order  to  avoid 
breakup  of  the  pulse,  the  pulse  is  subjected  to  frequency 
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discrimination  (outside  the  wiggler)  which  suppresses 
growth  of  sidebands.  The  frequency  discriminator 
effectively  acts  as  a  band  pass  filter  which  attenuates 
frequencies  above  and  below  the  pulse  frequency.  With 
stable  pulse  propagation,  a  stationary  state  is  reached 
when  the  energy  extraction  from  the  electron  beam 
balances  the  energy  losses  in  the  frequency  discriminator 
and  at  the  mirrors. 

The  electron  equations  of  motion  are  derived  from 
an  approximate  relativistic  Hamiltonian  in  which  the 
transverse  motion  is  taken  to  be  nonrelativistic  and 
transverse  spatial  variations  and  beam  self-fields  are 
neglected.  The  electron  current  driving  the  radiation 
field  is  obtained  by  summing  over  all  the  single  particle 
currents.  These  equations,  together  with  Maxwell's 
equations,  determine  the  temporal  and  spatial  evolution 
of  the  electromagnetic  pulse  in  the  wiggler. 

The  set  of  equations  governing  the  operation  of 
the  FEL  as  an  oscillator  have  been  formulated  and  they 
were  studied  both  analytically  and  by  the  use  of  a 
simulation  code. 

The  linear  eigenmode  analysis  leading  to  a 
calculation  of  the  linear  gain  as  a  function  of  the  FEL 
parameters  is  described  in  Appendix  B.  Therein  is 
discussed  the  limit  of  ultra-short  electron  beamlets  in 
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which  the  electron  pulse  length  l  is  very  short 

k 

compared  to  the  slippage  distance  k  L  /k  >  r  l  . 

w  w  s 

The  linear  gain  per  pass  has  been  calculated  for 
constant  parameter  wigglers,  variable  parameter  wigglers, 
and  a  combination  of  constant  and  variable  parameter 
wigglers.  The  analysis  included  the  effects  of  frequency 
discrimination  as  well  as  a  finite  energy  spread  (modeled 
by  a  Lorentzian  distribution  function). 

It  is  found  that  the  gain  depends  on  only 
two  combinations  of  FEL  parameters  and  graphs  of  the 
linear  gain  as  a  function  of  these  FEL  parameters  have 
been  plotted  and  the  areas  of  good  linear  gain  located. 

A  notable  feature  of  the  linear  gain  is  its  dependence 
on  the  time  delay  Atre^  of  the  electron  pulse  with 
respect  to  the  time  of  entry  of  the  EM  pulse  into  the 
the  wiggler.  The  linear  gain  is  zero  when  At  is 

zero,  reaches  a  maximum  when  Atrel  >  0  is  finite,  and 
finally  decreases  with  increasing  Atre^.  These  graphs 
have  been  used  to  select  appropriate  parameters  for 
input  into  the  simulation  code. 

A  one-dimensional  code  has  been  written  to 
simulate  operation  of  the  FEL  in  which  the  electro¬ 
magnetic  pulse  is  grown  from  noise  levels  all  the  way  to 
saturation.  The  details  of  the  code  are  described  in 
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Appendix  C.  This  code  enables  us  to  study  both  pulse 
formation  and  the  stability  of  the  saturated  state. 


Due  to  computer  time  limitations,  most  of  the 
simulations  have  been  done  using  5-function  electron 
beams.  The  results  of  the  simulations  are  presented  and 
discussed  in  Appendix  C. 

In  a  typical  run  in  which  the  EM  pulse  grows  from 
noise  levels  to  saturation,  the  early  phase  of  the 
evolution  is  characterised  by  growth  of  the  linear 
eigenmode  with  the  pulse  amplitude  peaked  at  the  back. 

As  the  evolution  enters  the  nonlinear  phase,  the  maximum 
in  the  pulse  amplitude  moves  to  the  front  which  then 
continues  to  grow  at  approximately  the  linear  rate. 

Growth  continues  until  the  peak  amplitude  at  the  front  of 
the  pulse  is  large  enough  to  begin  to  trap  the  beam 
electrons.  At  saturation,  without  frequency  discrimina¬ 
tion,  the  pulse  remains  peaked  at  the  front  and  never 
broadens  towards  the  back.  The  electrons  are  not  effec¬ 
tively  trapped  throughout  the  length  of  the  wiggler  and 
energy  extracted  from  the  electrons  is  only  a  small 
fraction  of  that  which  is  theoretically  possible. 

The  failure  of  the  EM  pulse  to  broaden  towards 
the  back  is  due  to  the  growth  of  unstable  sidebands  (of 
the  order  of  the  electron  bounce  frequency  in  the 
ponderomotive  potential  well),  which  leads  to  a  loss  of 
coherence  and  eventual  detrapping  of  the  electrons. 


When  frequency  discrimination  is  introduced  by 

passing  the  EM  pulse  through  a  bandpass  filter,  the  final 

phase  of  the  evolution  depends  critically  on  the  band 

width  w  of  the  filter.  If  is  small  so  that  the 

F  F 

growth  of  unstable  sidebands  are  suppressed,  the  pulse 
broadens  into  a  square  pulse  and  the  electrons  trapped  in 
the  ponderomotive  potential  well  remain  trapped  during 
its  passage  through  the  wiggler .  The  saturated  pulse 
propagates  coherently  and  electron  trapping  is  very 
efficient. 

The  presence  of  frequency  discrimination  by  a 

band  pass  filter  has  an  effect  (proportional  to  l/w_) 

F 

which  tends  to  retard  the  advance  of  the  EM  pulse 
relative  to  the  electron  beam,  that  is,  it  tends  to 
reduce  the  effective  EM  pulse  advance  produced  by  a 
finite  Atj.^  >0.  To  avoid  the  EM  pulse  running  either 
ahead  of  or  lagging  behind  the  electron  beamlets,  it  is 
indicated  that  a  choice  Atre^o>F  *  1  is  desireable. 

This  was  confirmed  by  the  numerical  simulation.  Thus  to 
obtain  finite  linear  gain,  cannot  be  altered 

independently  of 

A  small  value  of  is  desirable  to  suppress 

the  sidebands  and  guarantee  a  stable  saturated  state. 
However,  small  values  of  u>F  typically  implies  small 
linear  gain  since  At  ^  ■  l/wF  is  then  large. 
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In  order  to  insure  that  a  value  of  ujp  may  be 
found  which  allows  both  for  linear  gain  and  stable 
propagation,  it  has  been  deduced  that  the  following 
minimum  current  criterion  must  be  satisfied: 


<I> 


> 


1.42  x 10 


5 


k  2r  2 
w  o 


(1  -  r)  Ay 

r 


Amps 


where  <I>  is  the  beam  current  averaged  over  the 

slippage  distance,  r  is  the  mirror  reflection 

coefficient,  rQ  is  the  beam  radius,  the  wiggler  length 

L w  »  ksrg  is  set  e<3uai  to  the  diffraction  distance,  and 

e2B2/m2k2c-  *  1. 
w  w 

Linear  gain  can  be  enhanced  somewhat, 
particularly  for  large  Atrg^,  by  adding  a  constant 
parameter  section  to  the  front  of  the  variable  parameter 
wiggler.  The  constant  parameter  section  cannot  be  too 
long  since  electron  trapping  may  be  adversely  affected. 
The  linear  growth  phase  in  the  constant  parameter  section 
will  tend  to  saturate  at  an  amplitude  inversely  propor¬ 
tional  to  the  square  of  the  length  of  the  constant 
parameter  section,  and  this  amplitude  may  not  be  large 
enough  to  trap  the  beam  electrons  efficiently. 
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In  view  of  the  above  considerations,  we  have 
simulated  an  FEL  with  the  following  physical 
parameters 


Electron  Energy  50  MeV 

Electron  Beam  Current 

averaged  over  slippage  distance  50  Amps 

Wiggler  Wavelength  1.25  cm 

Number  of  Wiggler  Wavelengths  400 

Signal  Wavelength  \  »  2 it w  /c  lp 

s 

Wiggler  Field  Amplitude  7.0  kG 

Nominal  Trapped  Electron  Energy  Loss  20  MeV 
Constant  Parameter  Length  7.5  cm 


Stable  propagation  of  the  saturated  state  was  observed 

with  a  filter  band  width  of  0.1  .  An  effective 

s 

trapping  of  ~80%  of  the  beam  electrons  occurred, 
resulting  in  an  energy  transfer  efficiency  of  ~30%. 
Without  frequency  discrimination,  the  energy  transfer 
efficiency  is  only  0.8%. 
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MAN LEY -ROWE  RELATIONS  FOR  THE  THIN  LENS  GAIN  EXPANDER 


BY 

M.  ROSENBLUTH  AND  V.  WONG 


I.  INTRODUCTION 

In  this  note  we  argue  that  two-dimensional  free 
electron  lasers  such  as  gain  expanders  may  be  described  by 
an  action-angle  Hamiltonian  of  the  form 


K  =  ko  J  + 


10) 

ee  s 


(t  +  (1/B(|  -1)  z) 


n,m 


(J,y) e 


i  (nk  z  +  m<}>) 


w 


where  J  is  the  phase  space  area  for  betatron  oscillations,  <t> 
its  conjugate  angle,  y  the  electron  energy,  i.e.,  E  *  ymc2 
and  x  *  t  -  z,  its  conjugate  momentum,  z  is  the  inde¬ 
pendent  variable  and  e  is  proportional  to  the  laser 
field. 

In  Section  II,  we  demonstrate  a  series  of  canonical 
transformations  which  explicitly  and  exactly  reduces  the 

1 

"thin  lens”  expander  as  recently  proposed  by  Madey,  et  al. 
to  the  above  form.  The  details  are  algebraically  messy  but 
the  procedure  is  straightforward.  In  Section  III,  we  show 
that  this  Hamiltonian  implies,  to  order  e2,  the  Manley-Rowe 
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relations,  and  generalized  Madey  theorems.  These  may  impose 

severe  restrictions  on  gain  expander  operation.  We  also 

display  in  terms  of  the  h  's  the  equations  for  the 

n,m 

energy  transfer.  In  Section  IV,  we  evaluate  these  for  the 
thin  lens  system  at  low  excitation  amplitude.  The  casual 
reader  should  restrict  himself  to  Section  III. 


II.  REDUCTION  OP  THE  THIN  LENS  GAIN 
EXPANDER  TO  NORMAL  FORM 

We  use  as  a  Hamiltonian  for  the  thin  lens  model  the 

z  momentum,  and  represent  the  wiggler  by  a  vector  potential 

A  .  So  for  y  >:>  1  we  have 
z 


K 

1 


(p-A) 2+l 
2y 


xF  (z) 
1 


X  ^ 

—  F  (z) 
2  2 


(1) 


z  is  the  independent  variable.  Y  may  be  regarded  as  a 
canonical  variable  and  x  =  t-z  as  its  momentum. 

For  the  thin  lens,  which  consists  of  a  wiggler 
composed  of  pairs  of  very  thin  magnets  with  transverse 
gradients,  we  have: 


(Z) 


F  (Z) 
2 


a  +  a 

+ 


[6. +  8 


6  ( Z-nX  ) 

o 


6 (Z-nX  ) 

o 


(2) 


A  is  the  vector  potential  of  the  laser  A  =  eA^/c. 
Since  electrons  with  no  transverse  excitation  pass  through 
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the  magnets  at  a  fixed  x  value,  the  coefficients  in 
equation  (2)  are  determined  by  expansion  around  these  points. 
Analytic  solutions  are  known  both  for  the  average  (nominal) 
motion  xq(z,y);  pQ(z,Y)  and  for  the  betatron  deviations 
from  average  (when  the  laser  is  not  present) .  We  take 
advantage  of  this  to  perform  a  series  of  canonical 
transformations,  which  reduces  the  Hamiltonian  to  the 
standard  form. 

The  first  transformation  involves  a  generating 

function 


G  (p',x,t',y)  »  p'x  -  x  p'  +  xp  +  t'y 

1  0  0 


The  usual  canonical  transformation  equations  are 

x'  =  9G/9p' 

P  *  9G/9x 
Y'  =  9G/9T ' 

T  =  9G/9Y 
K2  *  +  9G/9z 
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so  that 


x '  =  x  -  x 


P'  =  P  -  PQ 

3x  3p 

T'  =  T  +  p'  — -  (x1  +  X  )  — — 

Zy  o  3y 

Y  '  =  Y 


(3) 


and  the  new  Hamiltonian  is 


K 

2 


..  2 


— —  +  -  x'2  F  +  -2- 

2y  2  2  2y 


Note  that  for  p'  =  x'  e  0,  (no  betatron  oscillations) 
we  find 


dx 

dz 


1 


Proper  choice  of  a  and  8  makes  the  system  isochronous, 

i.e.,  3/3y  (<3t  /dz)  «  0.  We  assume  this  hereafter.  Next  we 

o 

introduce  a  second  transformation 
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G2  =  p"x' 


+  T  "y 


Yz 


which  simplifies  the  Hamiltonian  to 


K 

3 


-  EA  t|i 

Y  2y 


(4) 


with 


T "  =  T 


Z  +  p  ’ 


dx 
_ _ o 

dY 


(5) 


The  next  transformation  reduces  the  Hamiltonian  to 
harmonic  oscillator  form.  Here  we  take  advantage  of  the  fact 
that  the  harmonic  oscillator  coordinates  P,  Q  must  be 
linear  functions  of  the  coordinates  p',  x'.  Hence  we  choose 
a  generating  function 

P2g  (z)  ,2 

G 3  =  - ^ ~ -  g2(z)  +  Px'  g3(z)  +  t'"y  (6) 


yielding 
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Q  =  P  g  +  x '  g 

1  3 

p '  =  g^  x 1  +  P  g  ^  .  (7) 

The  functions  g  (z)  may  now  be  determined  by 

1  f  2  #  3 

requiring  that  the  above  equations  are  satisfied  for  the 

desired  harmonic  oscillator  form,  i.e.,  P  =  -  sin  k.z, 

p 

Q  =  cos  kgZ  and  where  x',  p*  are  a  known  thin  lens 
solution  as  derived  previously: 


x'  =  3  f  1-*-  (  — )  nnl  cos  nk  A  +  —  (z-nA  ) 
L  J  Boy  o 


These  may  easily  be  checked  from  the  equations  of  motion. 
Here  a  is  a  normalization  coefficient,  and  we  are  at 

position  (n+1)  A  >  z  >  nA  . 

oo 

In  terms  of  the  magnet  parameters,  we  have  defined: 

(9a) 

with  8_>0,  and  the  betatron  frequency  satisfies 
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cos  k„  A 

S  o 


'8-M2 


By  substituting  (8)  into  (7)  and  laboriously  matching 
coefficients  of  cos  k  z,  sin  k  z,  we  find: 

P  P 


gi  =  A(z)/B(z) 


g2  =  C(z)/B(z) 


where 


g  =  1/B (z) 


—  —  - - 5—  )[l+(-)  nl  sin  k  z 

(sin  k  X  )  (1-n  )Y  )  l  j  B 

Bo  / 


j~  (l- (-)  nn)  sin  kg  (z-A  )  *  [l+(-)nn]  sin  > 

0  L  0  1) 


(sin  k  X  )  (1-n2  )Y 
B  0 


[i+(-)nn] 


cos  kgZ 


) 

+  —  |(l—  ( - ) nn)  cos  kg  (z-X^)  -  (l+(-)nn)  coskgzj| 


C  =  /- - : —  -y  - 5—  <  (l- (-) nn)  cos  k  (z-A  ) 

VXoSlnVo(1-n  1  80 


-  (l+(-)  nn)  cos  k  z 


27 


Here  z 


=  z  -  nAo  and  the  f's  are  hence  proper  periodic 
(period  2Aq)  functions  of  z. 

A  detailed  calculation  of  the  new  Hamiltonian  may  now 
be  made  using  the  equations  of  motion  and  after  vast  algebra 
we  find  as  we  knew  must  happen: 


K 


+  Q2) 


pA  A2 
—  +  — 

Y  2y 


A  final  transformation  to  action  angle  variables  is 
performed  involving  a  generating  function: 


G  =  T11"  y  +  J 


- 1 

sin 


JL_  + 

✓2J  /5j 


V  1  -  Q2/2J 
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III.  GENERAL  THEOREMS  FOR  STANDARD  HAMILTONIAN 


Thus  our  final  Hamiltonian  is 


K 


V“ 


2A  +  Al 

Y  2y 


(ID 


with  J  =  (P2  +  Q2)  /2. 

A  price  has  been  paid  for  this  simple  form.  The 
canonical  variables  are  J,  its  conjugate  angle  <f>,  y,  and 
a  transformed  time  timi  i  t.  The  added  complication  in 
the  time  variable  arises  from  the  implicit  dependence  of 
the  f's  since  the  timelike  variable 


A 

T 


3G 

3 

3Y 


(12) 


Note  that  y  has  survived  as  the  canonical  coordinate. 

Before  proceeding  further  it  is  worthwhile  to 
determine  the  physical  significance  of  J  =  (P2  +  Q2)/2.  By 
using  (7),  (9)  and  (10)  we  can  compute  J  as  a  bilinear  form 
in  the  p' ,  x'*s.  Thus  the  constant  J  surface  is  a  tilted 
ellipse  in  the  original  betatron  coordinates.  We  find  as 
expected  from  the  area  preserving  properties  of  the 
canonical  transformations  that  2ttJ  is  simply  the  area  of 
that  ellipse.  This  is  consistent  with  the  observation  that 
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the  Hamiltonian  (11)  would,  if  were  a  slowly  varying 

function  of  z,  have  J  as  an  adiabatic  invariant. 

Let  us  now  turn  our  attention  to  the  laser  field.  We 
may  write  the  Hamiltonian  (11)  in  the  form  (expanding  in 
laser  amplitude) 

H  +  H 

o  1 

V 

where  for  a  plane  wave  laser  propagating  in  the  z  direction 
with  vector  potential 

A  «  a  cos  oj  (t-z/c)  (13) 

s  s 


we  find 


H 

l 


p'+P  (Z)  10)  T 

0  s 

-  e 


+  complex  conjugate. 


(14) 


Here  e  *  as/2. 

The  quadratic  term  (A2)  is  irrelevant  to  our 
purposes,  t,  p',  and  x'  must  of  course  be  expressed  in 
terms  of  the  canonical  coordinates. 
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Using  the  redefinition  of  variables  we  have 
introduced  along  the  way  we  see  that  Hi  may  be  written  as: 


H  =  £  [p  (z) 

1  Y  L  0 


lCu 

s 

e 


+ 


+  /j"  a(z)  cos  <$>+B(z)j  x 

z  +  ip(z)  y/j  cos  (z)) 
pCz)  J  cos  (2<j>+A(z))  +  £(z)J 


Here  all  the  functions  of  z  may  also  depend  on  Y,  are 
periodic  with  period  2*q,  and  could  be  determined  from  the 
foregoing  arithmetic  with  some  effort.  However  in  order  to 
see  the  nature  of  the  solutions  it  is  necessary  only  to  note 
that  we  may  write  our  Hamiltonian  in  the  general  form 


K 


£hn.m(J'Y) 


ink  z+m<i> 
w 

e 


+  complex  conjugate.  (15) 


While  we  have  derived  this  only  for  the  special  thin  lens 
case  it  should  be  generally  applicable  for  different 
magnets.  Note  that  any  x  dependence  of  the  laser  field 
would  not  affect  the  form  of  this  Hamiltonian. 
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The  point  of  all  this  maneuvering  is  that  we  may 
write  the  Liouville  equation  for  the  distribution  function 
f  (J,  y,  t,  4>)  s 


3f 

dz 


dJ 


♦  H  [(H)  *] 


(16) 


0  . 


Expanding  in  powers  of  e  we  have  3K  /3t  =  3K  /3 *  =  0. 

o  o 

We  assume  that  f  is  independent  of  t  (optical  phase)  and 
<)>  (betatron  phase).  The  latter  assumption  (certainly  applying 
to  a  case  with  no  initial  excitation)  could  probably  be 
relaxed  for  our  purposes.  Then 


3f  3K  3f  3K  3f 

_ i_  _ l  _ o  _ _ o 

3z  3  <p  3j  3t  3y 

3K  3f  3K  3f 

_ o  _ ]_  _  _ o  _ i 

3  J  3$  3y  3t 


(17) 

=  0 


and  if  we  write 


ICO 


fi  “  e 


HD'1  Z) 


+  complex  conjugate, 


we  have 
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We  further  keep  only  terras  which  are  not  oscillating  rapidly 

on  the  scale  of  k  to  find: 

w 


-iA  z 
n,m 


A 

n,m 


U) 

s 


+  complex  conjugate 


where 


A 

n,m 


nk  +  w 
w  s 


(f-) 


+  m  k 


B 


a) 

s 


(19) 


(20) 


sin2  A  z/2 
_ m,n _ 

(A  /2  \ 2 

\  n,m  J 


3f 
_ _ o 

3  J 


■j 


s 
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We  see  immediately,  by  multiplying  by  J  and  integrating  by 
parts  over  dy  dJ,  that  the  mean  action  change  is  given  by: 


o> 


-  j  JfodJdl""["  h  -  “s  £]  |hn. 


m 


(21) 


sin2  A 


m 


z/2  ) 
-n _ l 


(V)!  ) 


and  the  energy  change  by: 


<5y)  =  - 


u >  e 

s 


fa  a  ‘ 

m  —  -  to 

!  h 

s  3y 

|  n,m 

sin2  A  z/2  I 

_ m,n  V  (22) 

<V)!  i 

In  the  usual  case  where  only  one  resonance  is  of 
importance  (which  is  practically  necessary  and  desirable  to 
get  cancellation  of  excitation  at  a  defintie  z  value)  the 
Manley-Rowe  relation  follows  directly: 
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We  have  already  remarked  that  2ttJ  is  just  the  phase  space 
area,  hence  can  not  be  changed  by  moving  from  one  magnet  to 
another.  Only  going  from  one  m  to  another  seems  to  offer 
the  possibility  of  excitation  cancellation  with  gain  but  we 
haven't  yet  been  able  to  implement  a  viable  scheme. 

The  m  =  0  resonance  does  not  produce  excitation  but 
since  the  whole  effect  depends  on  energy  derivatives, 
presumably  large  energy  spread  results.  It  should  also  be 
noted  that  only  m  =  1  is  operative  for  small  J,  hence 
only  h i  i  need  be  calculated  which  is  doable  in  spite  of 
the  complications  (see  Section  IV) . 

By  taking  the  second  moments  of  equation  (20)  it  is 
also  easily  possible  to  generalize  the  Madey  Theorem  relating 
gain  to  spread  without  making  the  resonant  approximation, 
i.e.,  keeping  all  m  values: 

(5J>  =  2  Tj  ((5j)2>  +  \  ^  (fJ6Y) 

(24) 

<iY>  .  i  x  <(4t)'>  .  i  x  <jJ6Y>  . 
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IV.  ENERGY  GAIN  FOR  THE  THIN  LENS  EXPANDER 

In  order  to  calculate  energy  gain  at  small  excitation, 
J,  coming  from  the  resonance  m  =  n  =  —  1,  it  is  necessary 
only  to  keep  the  terms  in  H  proportional  to  e-^  and  to 
Jj.  From  equations  (13)  and  (5)  with  e  =  a/2,  we  have: 

S 


H  = 

l 


io>  p 

s  o 


f  dx  dp 

|_p  17  17 J 


Recalling  that  P  =  -  /2j  sin  <p ,  Q  *  /2J  cos  <J>  and 
using  equations  (7)  and  (8),  we  find  readily  that  the  terms 
in  p',  x'  which  are  proportional  to  e_i<t>  are: 


p'  =  VJ/2  [c  +  i  dJ 
x'  =  >/j/2  [b  +  i  a] 

where  A,  B,  C  are  defined  in  equation  (10)  and 
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We  further  need  the  properties  of  the  nominal  isochronous 
trajectories 


dpo 

Po  d7 


Y 


dx 

o 

Po  dY 


1  ^ 

2  8+  Y2 


(-) 


n 


(27) 


(28) 


where  the  0's  specify  the  magnet  gradients  as  given  in 
equation  (2). 

Equations  (26)  -  (28)  may  now  be  substituted  into 
equation  (25).  We  must  still  perform  a  Fourier  Transform  in 
z.  Note  that  the  desired  resonant  Fourier  component  is: 
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1 
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i 


,  - 1 


X 

o 

/+  i  TT  z/X 
e  o 

-X 

o 


H  (z)  dz 


(29) 


and  that  the  integrand  consists  of  functions  of  z  =  z  -  nX^ 
times  constant  factors  and  factors  proportional  to  (-)n. 
Only  the  latter  contribute  to  (29)  .  It  is  thus  useful  to 
write  down  the  parts  of  equation  (26)  which  are  constant  and 
those  which  go  like  (-)n 


C.+  iD.  = 


ik„Z  I",  -ik  X  1  f 
-tie  6  j^l+e  B  oj  «/ 


c+  +  i  D+  =  elV  [e-Vo  .x]  J 


X  sin  k  X  (1-n2 ) 

o  Bo 


i  A _  +  B_  =  -  n/2  e 


ik0z 


(l  -  2z/Xo)  (] 


2z/X  )  (l  +  e‘lk6Ao 

o 


(l-e^Vo) 

)]f  v 


X  sink  X  (1-n2) 

Bo 


Here  the  subscript  -  denotes  the  terms  proportional  to  (-)n 
The  integrals  are  now  all  elementary  so  that  we  find 
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with 


X  = 


X  I. 


CO 


(*+Vo}  2n(1  +  COSk8Xo)  +  ~  Ao 


-2(1  +  cos  k_  A  ) 


2 

p  n 
0 


3  o  (tt  +  VA  ) 
S  o 


+  sin  k  A 

8  o 


n+i; 


-H 


[(sin 


k8  V  (1-n2)YXo 


We  may  simplify  this  somewhat  if  we  evaluate  at 
resonance  where  u.g  Aq  (Pq2  +  1)/2y2  *  n  +  kg  Xq  to  obtain 
finally: 


X  = 


2rj_ 


(1+p 


o>  (  "  + 


1  +  cos  kn  A  ) 

_ 6_  0/  2.  .  , 

(1-p  )  +  sink.  A 
o  8  o 


X 

3  o 


1  + 


B_  1 

S+n 


(sin 

k 

S 

A  )  (1-n2)  Y  A 
0  0 


\-H 


(30) 


and  finally  from  equation  (22) 
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with  ag  =  eAs/m  the  dimensionless  vector  potential  of  the 
laser  as  defined  in  equation  (13) ,  and  A  defined  in 
equation  (9) . 

The  term  proportional  to  8_/B+n  is  usually  the 
largest  in  equation  (30)  and  yields  the  same  result  as  Kroll, 
et  al.  The  m  =  +  1  resonance  result  is  obtained  by 
changing  the  sign  of  in  equation  (30) . 

Those  results  are  in  agreement  with  numerical 
simulations  by  John  Madey. 3 
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I .  INTRODUCTION 

In  the  free  electron  laser  (PEL) ,  a  relativistic 

electron  beam  is  propagated  through  a  transverse  periodic 

(wave  number  kw)  magnetic  field  wiggler  where  the  electrons 

acquire  transverse  oscillatory  motion  (amplitude  p  )  which 

enables  them  to  couple  to  the  transverse  electric  field  of  an 

electromagnetic  wave  (frequency  u>  and  wave  number  k  ).  For 

s  s 

low  density  beams  (where  collective  effects  may  be  neglected) 

the  coupling  may  be  viewed  in  terms  of  the  interaction  of 

electrons  with  a  ponderomotive  potential  well  (produced 

by  the  combined  wiggler  and  radiation  field)  moving  with 

velocity  v  =  w  /(k  +k  ) .  Electrons  with  velocity  close  to 
rss  s  w  s 

v  ,  or  equivalently  with  energy  close  to 
i  ©s 


are  in  resonance ,  and  net  energy  transfer  can  occur  between 

the  electrons  and  the  radiation  field.  Those  with  energy 

slightly  larger  (smaller)  than  y  lose  (gain)  energy  to 

res 

(from)  the  radiation  field. 

In  constant  parameter  wigglers  where  vres  is 
constant  throughout  the  wiggler,  the  energy  extraction  from 
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the  electrons  terminates  when  their  energy  is  reduced  to 
Yres*  This  puts  severe  limitations  on  the  gain  and 
efficiency  of  the  constant  parameter  FEL.  To  overcome  this 
limitation,  variable  parameter  wiggler s1  have  been  proposed 
in  which  electrons  are  trapped  in  the  ponderomotive  potential 
well  and  Yre„  gradually  reduced  by  adiabatic  variations  of 
the  wave  number  (k  )  and  amplitude  (p0  )  of  the  wiggler.  The 
energy  loss  of  the  trapped  electrons  translates  into  energy 
gain  of  the  radiation  field.  In  this  way,  the  energy 
extraction  in  a  variable  parameter  wiggler  can  be  made  to  be 
much  larger  than  that  in  a  constant  parameter  wiggler. 

In  principle,  a  DC  electric  field  may  be  applied  to  keep 
electrons  of  fixed  y  in  resonance  with  a  constant  parameter 
wiggler. 

While  a  variable  parameter  FEL  may  be  operated  as  an 
amplifier  by  passing  a  finite  amplitude  electromagnetic  pulse 
through  the  wiggler,  it  is  desirable  to  initiate  operation  of 
the  FEL  by  growing  the  pulse  from  noise  levels.  This 
necessitates  that  the  pulse  be  reflected  and  passed 
repeatedly  through  the  wiggler,  and  its  amplitude  increased 
on  each  pass.  However,  the  successful  operation  of  the  FEL 
in  this  mode  (as  an  oscillator)  will  require  that: 

1.  The  linear  gain  per  pass  be  sufficient  to 
overcome  reflection  losses  and  grow  the  electromagnetic 
pulse  to  a  large  amplitude  in  a  finite  number  of  passes; 
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2.  The  saturation  of  linear  gain  occurs  at 
amplitudes  large  enough  to  trap  a  significant  fraction  of  the 
beam  electrons  in  the  ponderoraotive  well; 

3.  The  finite  amplitude  saturated  pulse  propagates 
without  deterioration  due  to  the  onset  of  trapped  electron 
sideband  instabilities.1 

An  investigation  has  been  undertaken  to  determine  the 
extent  to  which  these  requirements  can  be  optimized  and  inte¬ 
grated  into  an  efficiently  functioning  variable  parameter 
PEL.  The  investigation  will  involve  not  only  theoretical 
analysis  of  the  PEL  equations,  but  also  numerical  simula¬ 
tions,  particularly  when  studying  the  nonlinear  saturation  of 
liner  gain,  electron  trapping,  and  the  effects  of  trapped 
electron  sideband  instabilities.  This  paper  is  primarily  an 
analytic  study  of  linear  gain.  In  a  subsequent  paper,  the 
numerical  simulations  will  be  described,  and  the  formation  and 
stable  propagation  of  the  finite  amplitude  pulse  discussed. 

Our  numerical  simulations  are  similar  to  those  in  the 

2 

pioneering  work  of  Colson  on  optical  pulse  formation, 
although  we  extend  his  techniques  in  order  to  treat  the  case 
of  interest  to  us,  the  high  extraction  variable  parameter 
wiggler . 

In  Section  II,  the  set  of  equations  governing  the 
operation  of  the  FEL  are  briefly  formulated.  A  fuller 
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discussion  may  be  found  in  Reference  1.  The  electron 
equations  of  motion  are  derived  from  an  approximate  rela¬ 
tivistic  Hamiltonian  in  which  the  transverse  motion  is  taken 
to  be  nonrelativistic,  and  transverse  spatial  variations  and 
beam  self-fields  are  neglected.  Thus  our  analysis  is  one¬ 
dimensional.  While  rough  estimates  show  that  two-dimensional 
effects  such  as  diffraction  and  radially  dependent 
ponderomotive  wells  are  not  crucial,  much  further  work  is 
required  in  this  connection. 

These  equations,  together  with  Maxwell's  equations, 
determine  the  temporal  and  spatial  evolution  of  the  electro¬ 
magnetic  pulse  in  the  wiggler.  In  order  to  suppress  the 
growth  of  sideband  instabilities  and  ensure  stable  propaga¬ 
tion  of  the  finite  amplitude  pulse,  the  formulation  includes 
frequency  discrimination  of  the  pulse  outside  the  wiggler, 
the  effect  of  which  is  to  attenuate  frequencies  above  and 
below  the  desired  pulse  frequency  w  .  As  we  will  see  in 
this  paper,  while  frequency  discrimination  is  necessary  to 
suppress  the  sideband  instability,  it  may  have  a  detrimental 
effect  on  linear  gain.  In  a  similar  fashion,  a  long  constant 
wiggler  section  will  be  found  to  be  desirable  for  linear 
gain,  but  to  lead  to  saturation  at  too  low  a  level  for 
trapping.  Thus,  careful  design  is  required  for  a  successful 
high  extraction  FEL. 
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The  FEL  equations  are  linearized  and  the  linear 
eigenmode  analysis  leading  to  a  calculation  of  the  linear 
gain  as  a  function  of  the  FEL  parameters  is  described  in 
Sections  III  and  IV. 

The  special  case  of  ultra-short  electron  beams  is 
considered  in  this  paper.  By  an  ultra-short  pulse,  which  we 
approximate  in  the  paper  by  6(z-vt),  we  mean  one  such  that 
the  difference  in  electron  and  photon  velocity  is  large 
enough  that  a  photon  is  able  to  slip  through  the  entire 
electron  beam  while  traversing  the  wiggler.  The  motivation 
(aside  from  the  intrinsic  interest  of  short  pulses)  for 
considering  this  case  are:  (a)  Analytic  tractability  and, 
more  important,  (b)  the  number  of  particle  orbits  which 
must  be  followed  in  numerical  simulation  is  obviously  far 
fewer  for  a  single  electron  beamlet  than  for  the  many 
beamlets  which  would  be  required  to  simulate  a  long  con¬ 
tinuous  pulse. 

The  basic  physics  of  electron-photon  interaction 
should  be  qualitatively  the  same  for  a  long  pulse  whose 
current  times  slippage  length  equals  the  number  of  electrons 
in  our  6  function  pulse.  The  linear  gain  for  long  electron 
pulses  will  be  considered  in  a  companion  paper. 

The  analysis  encompasses  constant  parameter  wigglers, 
variable  parameter  wigglers,  and  a  combination  of  constant 


and  variable  parameter  vigglers,  with  and  without  frequency 
discrimination. 

In  Section  V,  we  summarize  our  results  for  linear 
gain  and  present  a  semiquantitative  discussion  of  nonlinear 
propagation  and  sideband  suppression. 


II.  BASIC  EQUATIONS 


A.  Electron  Equation  of  Motion 


The  electron  equations  of  motion  can  be  derived  from 
1 


the  Hamiltonian 


H  =  H[x,  Px,  y,  Py,  (-£),  t,  z] 


c  Az  (1) 


where  the  electron  energy  (-£)  plays  the  role  of  the 
momentum  conjugate  to  t,  and  z  is  the  independent 
variable. 

4*  is  the  electrostatic  potential  of  an  accelerating 

9*o 

electric  field  -  — —  .  The  beam  density  is  assumed  small  so 

dz 

that  the  electrostatic  self-fields  of  the  beam  may  be 
neglected. 

The  vector  potential  A  is  the  sum  of  a  circularly 
polarized,  periodic  magnetic  field  wiggler  Aw: 
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a 

— w 


and  a  circularly  polarized  electromagnetic  plane  wave  A 

~  S 


=  a  (z,t)  x  cos(k  z  -  a)  t  +  C[z,tl) 
s  i  s  s 

-  v  sin(k  z  -  u>  t  +  C  [z  rt) )  I 
s  s  J 

where  aw,  kw,  ag,  and  C  are  slowly  varying  functions  of  z 
and  t.  As  has  been  shown  elsewhere  a  linearly  polarized 
wiggler  gives  only  slightly  modified  equations. 

The  assumption  that  A  is  independent  of  transverse 
coordinates  (x  and  y)  implies  that  the  canonical  perpen¬ 
dicular  momenta  Pv,  are  constants  of  the  motion.  We 

x  y 

will  assume  P  ■  0,  P  *  0,  an  assumption  corresponding  to 

x  y 

the  neglect  of  transverse  velocities  in  the  incoming  beam 
before  it  enters  the  wiggler. 
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If  we  introduce  as  new  variables  the  energy  parameter 


(?-  <%) 


Y  me2 


and  the  relative  phase 


k  dz  -t-  k  z  - 

W  3 


the  Hamiltonian  H  takes  the  form: 


H  =  H  (-  g,  t,  z) 


=  -  me  [y2  ”  U2  +  2a  a  cos  (ip  +  C)  ] 


-  mey  + 


|u2  -  2awas  cos  ( ip  + 


where  u  =  1  +  «w  +  V  aw  »  ag,  and  y2  »  u2  »  2  aw  as  . 

The  electron  equations  of  motion  in  terms  of  Y  and 
<Pt  derivable  from  H{-  g,  t,  z) ,  are  then  given  by: 
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lx 

dz 


3H 

3t 


d<b 

e  o 
me2  d2 


k  a  a 
s  w  s 


sin(^  +  C) 


(3) 


k 

w 


+  k 

s 


+  0) 

s 


3H 

3? 


'  k 


w 


(4) 


where  <*is  -  kg  c#  and  terns  of  order  <<  i  _i_  _ E  , 

,  .  2  ws  3t  '  ws  3t 

•‘avas/tJ  <:<  1/  have  been  deleted. 

The  relative  phase  *  varies  slowly  for  electrons 
with  values  of  Y  close  to  Yr  (the  resonant  energy) 
defined  by: 


These  electrons  interact  strongly  with  the  radiation  field 
and  will  exchange  net  energy  with  it.  Electrons  with  values 
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of  Y  far  from  Yr 
field. 


have  negligible  effects  on  the  radiation 


It  is  therefore  convenient  to  describe  the  electron 
motion  in  terms  of  the  variables  Y  *  y  -  Yr  and  <P 

In  the  limit  of  Y  K<  Y r  ,  the  electron  equations  of 
motion  may  finally  be  approximated  by 


d£ 

dz 


dz  me2  dz 


k  a  a 
s  w  s 


sin(ip  +  S) 


(5) 


dj,  _  2  kw  l 

dz  y 

r 


(6) 


and  are  derivable  from  the  approximate  Hamiltonian 
(  Yf  <l>,  z)  : 


(  <*Y  d<p  ) 

+  )  r  +  e  o  ( 

|  dz  me2  dz  j  ^ 


k  a  a 
s  w  s 


cosfip  +  ;) 


(7) 


0 


d\ 

We  define  a  constant  parameter  wiggler  by  — —  = 

d<j>_  dz 

and 


d<j> 

o 

dz 


=  0. 

In  a  variable  parameter  wiggler,  Y  is  allowed  to 
dYr 

vary  spatially,  ?  0  .  The  presence  of  an  accelerating 

electric  field  d^^/dz  affects  the  electron  phase  space  (y»<M 

trajectories  in  a  manner  similar  to  spatial  variations  of 

Y  ,  as  may  be  seen  from  Equation  (7).  In  subsequent 

dY_ 

sections,  it  will  be  assumed  that  Yr  is  constant,  =  0  . 
A  constant  or  variable  parameter  wiggler  is  then  modeled  by 


d<p 


B,  Electromagnetic  Field  Equations 

The  transverse  electron  beam  current  density  Jx 
determines  the  time  and  spatial  evolution  of  Ag  through 
Maxwell's  equation: 


4  TT 

c 


J 

-i 


Substitution  of  the  assumed  form  of  Ag  yield: 
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(3) 


3a  3a 

_ s  _ s 

3t  C  3z 


2TTe 

k  me2  2 


a 

s 


2ire 

k  me2  1 
s 


(9) 


where 


J. 

3 


3  =  1,2 


~  /v 

e,  =  x  cos(k  z  -  w  t  +  S) 

- 1  s  s 

-  y  sin  (k  z  -  u>  t  +  C) 

s  s 

A  A  § 

e_  =  -  x  sin(k  z  -  u  t  +  C) 

2  S  S 

-  y  cos  (k  z  -  u>  t  +  C ) 

s  s 

The  assumption  that  ag  and  £  vary  slowly  in  space 
and  time  has  again  been  invoked  to  delete  higher  order  terms. 
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T 


If  )  is  the  distribution  function  of  the 

o  o  o 

electron  beam  on  entry  at  z  *  0  into  the  wiggler  at  time 
t  =  -  — ,  the  beam  density  n(z,t)  and  transverse  beam 

O  Uig 

current  density  Jj(z,t)  are: 


n(z,t) 


)  5  [t  -  t  (y  ip  ,  z)  ]  dy  dip 
o  o  o  o  o 


(10) 


J  .  (z  ,t) 
3 


F  (y  . 
o  o 


e  6  It  -  t(Y  .  P  ,  z)  ]  dY  dip 
— ]  o  o  o  o 


(11) 


where  the  transverse  electron  velocity  vx  is  (Px  *  0): 


v 

— i 


c 

—  a 
y  -w 


Si 


c 

Y_ 


a  (z) 

V) 


cos 


>(v 


z)  +?(z,  t (Y, 


z)  1 


V 

— 1 


—2 


% 


—  a  sinl^  +  C} 

Y  w  '  ' 

r 
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*<V  '^o'  z )  and  *(Yof  tyQ,  2)  are  the  phase  space 
coordinates  of  the  electron  trajectory  with  initial  values 
and  at  z  =  0.  t(Y0,  2)  is  the  time  of  arrival 

of  the  electron  at  z  and  is  given  approximately  by 


t(y  »  4*  *  z) 
o  o 


£ 

V 


where 


Let  F0(Y0,  Vq)  *  where  f(YQ) 

determines  the  energy  distribution  and  J  dy  f(y  )  =  l. 

o  o 

Then,  from  Equation  (10): 


n(z,t)  =  H  -  tj 


and  from  Equation  (11)  : 
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v*-*1  «(M 

JdY0  f(Yo)  cqs|*(yo.  *o,  zl  +t(z,t)| 

*o  * 

»hen  t  changes  by  At  .  ^  or  z  by  Az  =  2SE  .  the 

argument  *0  of  V { :<o ,  Yo ,  z)  changes  by  2t.  Assuming 
that  the  beam  density  H(§  -  t)  series  sidy  over  these 
email  changes  in  t  or  z,  it  is  appropriate  to  average 
cos  |*t  4(  over  Thus 

Jl(2't>  *  f  V*>  H(M 

r 


_1_ 
2  7T 


/• 


dip 

o 


f  ( Y  ) 
o 


cos 


{* +  c} 


eca  (z) 
w 


Similarly 
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J2(z,t)  =  - 


eca  (z) 
w 


H  (v‘  fc)  <Csin  + 


The  field  equations  for  as  and  £  are  then  given 
by 


c 


3a 
_ s 

3z 


=  - - —  H  /sin  { ip  +  c|> 

k  cmy  x  1  ' 

s  r 


a 


s 


+  c 


2ire2a 

W  / 

=  -  H  Ceos 

k  cmy  N 
s  r 


{*  +  «(> 


(12) 


(13) 


C.  FEL  Equations 

In  the  FEL  configuration  which  will  be  analyzed,  a 
combination  of  constant  and  variable  parameter  magnetic  field 
wigglers  is  positioned  between  mirrors  so  that  an  electro¬ 
magnetic  pulse  can  be  reflected  repeatedly  and  multiple 
passes  of  the  pulse  through  the  wiggler  achieved  (see  Fig.  1) . 
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Relativistic  electron  beamlets  of  finite  length  are  injected 
into  the  wiggler  at  periodic  intervals  so  that  on  each 
forward  pass  of  the  pulse  through  the  wiggler,  there  is 
overlap  of  the  pulse  and  beamlet.  The  pulse  grows  on  each 
forward  pass.  Eventually,  the  amplitude  of  the  pulse  becomes 
large  enough  to  trap  a  fraction  of  the  beam  electrons.  At 
this  stage,  the  pulse  is  susceptible  to  detrimental  spatial 
distortions  due  to  the  onset  of  trapped  electron  sideband 
instabilities.  To  avoid  breakup  of  the  pulse,  it  will  be 
necessary  to  suppress  the  growth  of  the  sideband  instabil¬ 
ities.  This  may  be  accomplished  by  passing  the  pulse  through 
a  frequency  discriminator  which  effectively  acts  as  a  band 
pass  filter,  attenuating  frequencies  above  and  below  the 
desired  pulse  frequency.  With  stable  pulse  propagation,  a 
stationary  state  is  reached  when  the  energy  extraction  from 
the  electron  beam  is  balanced  by  the  energy  losses  in  the 
frequency  discriminator  and  at  the  mirrors. 

To  facilitate  formulation  of  the  equations  which 
describe  the  operation  of  the  FEL,  it  is  advantageous  to 
transform  the  variables  z  and  t  to  the  variables 
u  and  v : 
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where  L  is  the  length  of  the  wiggler. 

The  line  z  »  0  in  z  -  t  space  maps  into  the  line 
v  =  -  u  in  u-v  space  and  z  =  L  into  v  =  -  u  +  1  (see 

Fig.  2) . 

In  u-v  space,  the  beam  electrons,  entering  the 

wiggler  at  z  *  0  and  exiting  at  z  *  L,  move  on  lines  of 

constant  u.  The  photons  of  the  electromagnetic  pulse, 

propagating  in  the  beam  direction,  move  on  lines  of  constant 

v.  During  the  n  pass  of  the  pulse  through  the  wiggler, 

the  pulse  interacts  with  the  electrons  when  the  values  of  u 

and  v  lie  in  the  area  defined  by  >  u  >  u^  , 

l_u>v>_u  where  the  difference  u  -  u.  =  u  is 

in  in  o 

determined  by  the  length  of  the  injected  electron  beam. 

In  terms  of  the  u-v  variables,  the  coupled  set  of 
equations,  (5),  (6),  (12),  and  (13),  which  describe  the 
evolution  of  the  electromagnetic  pulse  passing  through  the 
wiggler  together  with  the  electron  beam,  can  be  rewritten  as 
follows: 
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Figure  2.  Electron  and  Photon  Trajectories 
in  u  -  v  plane 


65 


T 


where 


dip 

3v 


Y 


(14) 


ii  _ 

av  “ 
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.  ★ 
ia  e 


■i'i' 


(15) 


3a 

3u 


i  n  h(u) 


(16) 


a*  is  the  complex  conjugate  of  a,  and 


2  k  L  k  La 

.  W  S  W 

a  -  -  -  a 

s  y  Y  s 

r  r 


2  k  L 
~  w 

Y  =  — -  Y 


2  k  L  3<b 
w  e  o 


Yr  me2  3v 
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with  Ayr  the  change  in  y  which  would  be  experienced  by  an 
electron  freely  accelerating  under  the  potential  .  h(u) 

is  a  form  factor  which  determines  the  beam  density  profile 
and  is  defined  by: 


n(z,t)  =  h|^  -  t  j 


h  (u) 


(17) 


C  2n 

where  >  h(u)  d  u  *  1  and  N T  is  the  total  numer  of 

Uln  ,  , 
electrons  in  the  beam  per  unit  area,  NT=  jdz  n(z,  t)  . 

The  dimensionless  parameter  n  is  given  by: 


8 ire 2  N_  k  L2  a2 
_ T  w  w 

y  me  V  y  2 


(18) 
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Thus,  for  the  nth  pass  of  the  electromagnetic  pulse 
through  the  wiggler,  the  electron  phase  space  (Y,  ) 

trajectories  are  determined  by  the  solutions  of  equations 
(14)  and  (15),  with  initial  conditions  Y  *  YQ,  <P  =  at 
v  *  -  u.  The  phase  of  electrons  entering  the  wiggler  is 
uncorrelated  with  that  of  the  electromagnetic  pulse.  Thus,  *PQ 
is  distributed  uniformly  between  0  and  2tt.  The  average 
over  i>0  is  indicated  by  <  >  .  The  pulse  amplitude  after 

interaction  with  the  beam  electrons  is  given  by: 

an(v,o)  =  an(v,  -Uq)  +  Aan(v)  (ly) 


where 


A  a  (v)  =  i  n 


o 

du  h(u)  <e-i<,(v'ul>  g(v,u) 
-u 


g(v,u) 


1 


o 


1  -  u  >  V  >  -u 


o 


>  u  >  -u 

o 


otherwise 


(20) 
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In  deriving  equation  (20),  the  u,  v  coordinates 

were  transformed  so  that  the  interaction  of  pulse  and 

th 

electrons  during  the  n  pass  occur  for  u  and  v  in  the 
range  0>u>-uQ,  l-u>v>-u.  It  will  be  assumed 
that  there  is  no  overlap  of  the  electromagnetic  pulse  from 
pass  to  pass,  and  consequently,  this  transformation  can  be 
carried  out  in  deriving  corresponding  equations  for  the 
(n  +  1) th  pass,  etc. 

The  electromagnetic  pulse  is  reflected  backwards  and 
then  forwards  for  the  (n  +  1)  pass  through  the  wiggler. 

n+i 

Without  optical  filtering,  the  pulse  amplitude  a  (v,  -u  ) 
before  the  (n  +  l)th  pass  may  be  related  to  the  pulse 
amplitude  an(  v,  0)  after  the  nth  pass  by 


^  ^  v  a  n .  . 

a  (v  -  B,  -  u  )  =  r  a  (v,o) 
o 


UQ)  +  A  an(v)J 


(21) 


where  r  <  1  accounts  for  the  reduction  in  amplitude  due  to 
energy  losses  on  reflection  at  the  mirrors,  and  6  repre¬ 
sents  the  tunable  pass  to  pass  shift  in  the  position  of  the 
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pulse  relative  to  the  front  of  the  electron  beam  at  the 
moment  of  entry  into  the  wiggler.  As  we  will  see  the  linear 
growth  rate  depends  sensitively  on  6  which  is  physically 
determined  by  the  distance  between  mirrors  and  the  interpulse 
spacing  of  the  e-beam. 

If  identical  electron  beams  are  injected  into  the 
wiggler  on  each  pass  with  fixed  periodicity,  the  pass  to  pass 
change  in  pulse  amplitude  is  determined  by  equations  (19)  and 
(21),  independent  of  the  value  of  n. 

To  suppress  the  growth  of  sideband  instabilities  when 
the  pulse  amplitude  is  large  enough  to  trap  electrons,  the 
pulse  can  be  filtered  to  attenuate  sideband  frequencies  above 
and  below  the  desired  pulse  frequency  ws  .  This  filtering 
may  be  accomplished  by  a  band  pass  filter,  modeled  by  the 
equation: 


^  /sn 

—  a  (v,o)  +  v  a  (v,o)  =  v  a  (v,o) 

3v 


where  a  is  the  filtered  pulse  amplitude,  and  the 
frequency  half  width  Au)  of  the  band  pass  filter  is 
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VC 


~n+l  -n 

Since  a  (v  -8,  -  uq)  =  ra  (v,0)  ,  the  pass  to 
pass  change  in  pulse  amplitude  after  interaction  with  the 
beam  electrons,  filtering,  and  reflection  at  the  mirrors,  is 
given  by 


1 

v 


3v 


»n  +  1 
a 


(v  -  6, 


u  ) 
o 


+ 


~n  +  1 
a 


(v  -  8, 


u  ) 
o 


(23) 


*  r  £an(v,  -  uq)  +  A  an  (v)j 


Equations  (14),  (15),  (21),  and  (23)  govern  the 
operation  of  the  FEL. 
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III.  LINEAR  THEORY— WITHOUT  FREQUENCY  DISCRIMINATION 


For  small  pulse  amplitudes,  the  FEL  equations  are 
linearized  to  obtain  the  linear  eigenmode  equations.  In  the 
limit  of  short  electron  beams  (uq  <<  1)  where 


1? 


equations  (14),  (15),  and  (21)  yield: 


V 

(v,u)^>  =  j  J"  dv'  (v  -  v' 


)  I (v  -  v' ) 


u 

ftn,  '  ,  n  f  ,  •  /  -ii|»  (v*  ,u'  )\ 

a  (V  ,-uo)+x—  J  du  (e  / 


-u 


.  r  ,2 

X  -  (V  2 


1  >  v  >  0 


(24) 
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n  +  1 

a  (v  -  8,  -  u  )  -  r  a  (v,  -  u  ) 


Tt*  ) 


n  +  1 

A 

a  (v,  -  u 


n 

-  r  a  (v,  -  u  > 
o 


irtl 

■'-u 


-2-f  du  <e'i*(V-  "’> 


where 


I  (v  -  v’ ) 


f(Y  ) 
o 


i  Y  (v*  -  v) 
e  o 


1  >  v  >  o 

otherwise 

(25) 


and  the  relative  pass  to  pass  shift  8  of  the  pulse  is 
assumed  to  be  small  permitting  the  indicated  Taylor  series 
expansion.  r  is  condidered  to  be  independent  of  u  and 
v. 

Let 
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where 


where 


~n  +  1 
a 


(v. 


u  ) 
o 


6 

e 


'n 


(v. 


u  ) 
o 


*  (1  +  6)  an  (v,  -  u  )  , 

o 

|6|  «  1 


6  is  the  linear  gain  and  phase  shift  per  pass. 
Equation  (25)  may  therefore  be  approximated  by: 


(1  -  r  +  6) 


in  -  6  —  a 
3v 


n 


o 


in 

u 

o 


(  du  <e-i*(v'ul> 
-u 

o 


1  >  V  >  o 


(26) 


r  is  assumed  to  be  close  to  unity. 

The  boundary  condition  on  a  at  v  *  1  is 


an  (v  =  1 ,  -uq)  =  0 
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since  photons  with  v  >  1  have  never  interacted  with  the 
driving  electron  beam. 

Equations  (24)  and  (26)  are  valid  only  in  the 
interval  1  >  v  >0.  However,  the  upper  limit  may  be 
extended  to  infinity  since  the  solution  of  these  equations  in 
the  interval  1  >  v  >  0  are  unaffected  by  the  solutions  in 
the  interval  v  >  1.  With  this  extension  of  the  interval  to 
°°  >  v  >  0 ,  the  solution  of  these  equations  may  be  obtained 
by  Laplace  transformation. 

Define  the  Laplace  transforms 

1“  ;  r  v2 

dve'pV  <e-i4"v’u)>  e"'2 
o 

f  i.  y2 

A(p)  =  I  dve  pVS(v,  -uQ)  e  2 
o 


00 


O 


The  superscript  on  a  will  be  deleted  except  when  it  is 
necessary  to  differentiate  variables  referring  to  different 
passes  through  the  wiggler. 


75 


Equations  (24)  and  (26)  are  Laplace  transformed 


0  (u) 
P 


1 

2 


du  '  0  (u  ') 
P 


(27) 


iTB  A 
3p 


+  Bp  A  -  (1  +  6-r) 


A 


±IL 

u 


o 

1 


du  0  (u) 
P 


=  S  &  (o,  -u  )  =  a 
o  o 


(28) 
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1 


Equation  (27)  is  solved  to  obtain 


0p(u) 


,  31 

A  _ E 

2  3p 


31 

-i  n  _ p 

2  3p 


(.29) 


and  substitution  of  0p(u)  in  equation  (28)  yields: 


ire 


3 

—  A 
3p 


+  6pA-  (1  +  6  -  r) 


A 


+  A 


‘  .  31 

-i  n  p 

2  3p 

e 


1 


a 

o 


(30) 


The  function  I  is  determined  by  the  initial  energy 
P 

distribution  f (y  )  of  the  beam  electrons.  If 

o 
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simple  when  f(YQ)  is  Lorentzian,  and  consequently,  the 
Lorentzian  form  of  I  is  used  in  the  subsequent  analysis 
unless  otherwise  stated. 

Equation  (30)  is  solved  to  obtain  the  solution  A  (p) 
which  vanished  at  00  in  the  right  half  p-plane.  The  pulse 
amplitude  a(v,  -  uQ)  is  then  evaluated  from  A(p)  by  the 
inverse  Laplace  transform 


a  (v,  -  u  ) 
o 


-i  Tv2 
2 

2-rri 


,  pv 

dpe  A  (p)  ,  1  >  v  >  o 


(32) 


where  C  denotes  the  contour  of  the  inverse  Laplace 
transform. 

The  eigenvalue  equation  is  obtained  from  the  boundary 
condition 


/\ 

a  ( v  =  1 ,  -  u  )  =0 
o 


which  in  terms  of  A(p)  may  be  written 
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r 


e  2  j  dpep  A (p)  =  0 
c 

A.  Constant  Parameter  Wiggler 

In  a  constant  parameter  FEL,  r  is  zero.  Thus 
equations  (30)  and  (33): 


where 


(1  -  r  +  6)  +  3  ( 7.  ,  +  i  y  ) 
tn 


6 


1 


P  +  Yth+iY 


(33) 


,  from 


(34) 


(3  5) 


(36) 


80 


and  C.  is  obtained  from  C  through  the  transformation 
from  p  to  5- 

If  (n82)  V3  <<  l,  and  the  exponential  is 
approximated  by 


<fYs 


JL 


-l  * 


-i 


/  nS2  \V3  l 

12/  C2 


(which  corresponds  to  neglecting  the  second  term  within  the 
square  bracket  of  equation  (24) J ,  equation  (34)  simplifies  to 


j 


L 


(3£.  -  2 5 L ) 


0 


(37) 


where 


are  the  roots  of  the  cubic  equation 


(38) 
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82 


(p  +  Y  fch  +  iY) 

g  =  rr/2 


and  to  rewrite  equation  (30)  as  follows: 


— —  -  iqA  +  iAA 
3q  H 


sr"1/2 


in 

2Tq2 


(40) 


where 


A 


(1  -  r  +  5)  +3  (y  +  iy) 

_ tn _ 

sr  V2 


The  solution  of  equation  (40)  which  vanishes  when  p 
tends  to  infinity  in  the  right  half  p-plane  is 
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oo  ( i— i }  i 

ia0  Sl(q)  [  ,  .  ~Sl(q  ] 
rrw  e  \  dq  e 


where 


SL(q) 


i  A<3+^17T 


£cxj(l-i) 


in 

2fq  2 


The  eigenvalue  equation  is  therefore  given  by 


dq  e 


r1/2  q  +  (q)r(L  ^  -S1(q') 
dq  e 


In  the  limit  of  an  infinitely  long  wiggler  where 
r  -  the  integral  in  equation  (41)  may  be  evaluated  by 
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deforming  the  contour  Cg  towards  the  Im  q  axis.  The 
contour  approaches  a  saddle  point  at  <1  ~  A  +  i  F1/2  f  an<j 

if  Re  A  >  0/  the  contribution  of  this  saddle  point  is 

AF  V2 

exponentially  ~  e  large.  This  implies  that 

equation  (41)  can  be  satisfied  only  if 


00  ( 1“  i ) 

f  -si(q) 

I  dq  e  =  0 

-oo(l-i) 


(42) 


where  the  contour  passes  above  the  origin.  Equation  (42) 
determines  the  eigenvalue  $  in  the  limit  AT1/2  -*•  °°  . 

If  the  exponential  in  s1(q)  is  again  approximated 
by 


e 


in 

2  Fq2 


1 


in 

2Tq2 


+ 


equation  (42)  reduces  to: 


°2  ( 1— i) 

1  <*q  e 

-°°( 1-i) 


iq  .  ,  a 

-  —? —  +  i  A  q  -  — 


=  0 


(43) 


t 


where 


a  2sr3/2 

For  a  -*■  0/  equation  (43)  may  be  approximated  by 


Equation  (43)  determines  A  as  a  function  of  a  and 
it  is  readily  verified  that  A(a)  is  a  solution  of  the 
differential  equation 


(44) 
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For  a  >>  1#  a  solution  of  equation  (44)  is  given 


approximately  by 


i  7T 


+  C. 


a"1/9 


+  ii 


iff) 


1/3 


-i  TT 

e  ®  a-1/3 


(45) 


where  CI  is  a  constant  to  be  determined. 

Equation  (43)  is  solved  numerically,  and  in  Figure  4, 
the  variable  Re[(l  -  r  +  5  +  3  Yth) (2/nB2 )  ] 

is  plotted  as  a  function  of  2Sr3/2/n  =  1/a.  For  small  1/a, 
the  curve  follows  equation  (45)  with  —  1.61  -  i  1.94 
and  decreases  to  zero  at  ^  =  1.42. 

C.  Combination  of  Constant  and  Variable  Parameter  Wiqqler 

In  the  case  of  a  magnetic  wiggler  with  a  constant 
parameter  section  at  v^  >  v  >  0  and  a  variable  parameter 
section  at  1  >  v  >  v  ,  the  eigenraode  equation  may  be 
obtained  by  an  extension  of  the  analysis  previously 
described. 

Thus,  in  the  constant  parameter  section  >  v  >  0, 
the  pulse  amplitude  is  (equation  (32)): 
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a (v,  -  u  ) 
o 


a  ( 0 


-  u 

o 


) 


j  =  1,2,3 


(q j  -  e)  r1/2  v 


( 3q j  -  2A) 


v1  >  v  >  0 

(46) 


where 


■( 


vth  +  1  v 


and  q j  are  the  roots  of  the  cubic  q3  -  q2 


In  the  variable  parameter  section  1 

-i  iMv,u) 

relevant  equation  for  <  e 


> 


is : 


A  +  i  a  « 
>  v  >  v1# 


0. 

the 


<e-i4»(v,u)> 


ir 


(V  -  vL) 


■if 


dv'  (v^  -  v  ') 


I(v-v’)  a  (v  '  ,u) 


r1 

;  (V-V^  ^  dv' 


I  ( v  -  v  )  a  ( v  ,  u) 


(47) 


v 


i  r 


(v'  -  v  )  2 


+  —  j  dv'(v-v')  I(v-v')  a(v',u)  e 
v. 
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The  corresponding  equation  for  A  is: 


—  -  iq  A  +  i.\  A 
3q 


where 


M-ln/2rq‘-1)  a 

Jisi  i  ,in/2rq2  ,WW  VV) 

Br  \e  -1/  \~  +  ~lp~/ 


ia (Vj 

r 


(48) 


bi  (Vj) 


r  f  a  „  “(Yth+i>r)  (v,-v'i 

-  —  I  dv  (v  -v  )  a  (v  , -u  )  e  cn  1 

2  /  1  o 

J  I A 


(42a) 


>.  rVi 


b  (v  ) 
2  1 


•t/ 


.  •  -  ,  •  .  -Tyth+iy)  (v  -  v' ) 

dv  a  (v  ,-u  )  e  tn  i 

o 


(49b) 


is: 


The  eigenvalue  equation  analogous  to  equation  (43) 


iMUftiirt 
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00  ( 1  — i) 


/ 


dq  e 
-°°  (1-i) 


-iq2/2  +  iAq-ot/q 


b, (v, ) 


a (v  , -u  )  -  2  i  a 


i  i  ' 


W  ) 
+  - - — 


i  o 


=  0 


(5?) 


This  equation  determines  A  as  a  function  of  3  and 
F  ”v  .  Equation  (30)  is  solved  numerically,  and  in  Figure  4  the 
|l-r+<$+8  7th |(2/nB2  is  plotted  as  a  function 


variable  Re 

3/2 

of  2ST  /n  »  1/a  for  different  values  of  : 


Figure  4.  Constant  and  Variable  Parameter  Wiggler  without 

3/ 

Frequency  Discrimination  -  Re(i<5  )  vs  2BT  2/n 


IV.  LINEAR  THEORY  -  WITH  FREQUENCY  DISCRIMINATION 


With  frequency  discrimination,  the  pulse  amplitude 
equation  (eg.  23)  is 


(l  +  6  +  ..W-~+ll  (a  -  8  Is- 
\  )  v  3v  I  3v 


02  32  -  ^ 

+  -r-  r — r  a  + 
2  3vz 


(51) 


=  r  |  a  +  Aa 


where  6  is  defined  in  Section  III  and  B  <<  1. 

The  analysis  will  be  restricted  to  the  case  of  v 

P 

This  choice  is  motivated  by  the  consideration  that 
Equation  (22)  indicates  that  the  filtered  pulse  is  effec¬ 
tively  an  average  of  the  preceding  Av«-f  i.e.,  it 

v 

corresponds  to  an  effective  retardation  of  the  light  pulse  by 
^  relative  to  the  electron  beam,  while  3  represents  a 
forward  shift  of  the  light  pulse.  Unless  the  two  are  nearly 
equal,  the  electron  pulse  and  photon  will  separate  after  many 
passes.  Numerical  work  has  verified  that  this  choice  is 
approximately  optimal. 
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Equation  (51)  then  simplifies  to 


B2  3 2 


(1  -  r  +  6)  a  -  — 


2  3vJ 


o 

j  du  <e-ilHv'u)> 


1  >  v  >  0 


When  v  <  0  or  v  >  1, 


(1  -  r  +  <5)  a  - 


ii  i!i 

2  3v2 


The  solution  satisfying  the  boundary  conditions  a  -*•  0,  v  ±« 


{2  (l-r+6) Kv 
!  B 


v  <  0 


(2 (l-r+6) }  v 
B 


v  >  1 


94 


where 


Re(l  -  r  +  6)  >  0.  At  v  ■  0  and  v  ■  1,  a  and 
3a/ 3v  are  continuous. 

The  procedure  of  Section  III  may  now  be  followed  to 
obtain  the  eigenvalue  equation.  Thus,  Equation  (52)  is 
Laplace  transformed,  and  the  equation  for  A(p)  in  the  limit 
of  short  electron  beams  is: 


p2-ir  +  2ir  T-p-r2„2 

‘  3p  ^  3p2 


.  2 (l-r+6)  . 

+  - p -  A 


JT  A  j  e  2  3p 


3  a  a 

+  Pa 


(2Cl-r+6)} 
iP  +  6 


a  (0 ,-u  ) 
0 


Equation  (53)  is  solved  to  obtain  the  solution  A(p) 
which  vanishes  in  the  right  half  p-plane.  The  pulse  ampli¬ 
tude  a(v,-u0)  is  then  evaluated  from  A(p)  by  the  inverse 
Laplace  transform 
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If 


(nS2)1* 


<<  1  and  the  exponential  is  approximated  by 


Equation  (54)  simplifies  to 


a(v,-u  ) 
0 


a(°*-uo)  ^ 

j=l , 2 , 3 , 4 


(n/32)  (C.-5,)v 


'jVMFrM  -«.] 


(55) 


where  are  the  roots  of  the  quartic 

52(5-63)2  -  5262  +  i  -  0  .  (56) 

The  continuity  of  a  and  3a/3v  at  v  ■  1  is 
satisfied  if 


(VS’)'  5  ,  HV 

e  ’yvw 


_ 31  *3 

fj-*.)  («!•*■)  ‘  S‘. 


(57) 


Equation  (57)  is  the  eigenvalue  equation 


If  5=0  (no  thermal  spread),  the  roots  of  the 
quartic  are 


As  in  Section  III,  it  may  be  conjectured  that  in  the  limit 
(n/B  )  -*•■»,  Equation  (57)  can  be  satisfied  if  two  roots 

(£  and  £  )  of  the  quartic  have  equal  positive  real  parts. 

1  2 

The  requirement  for  this  to  occur  is 

5  *  2e11T/4  +  a  -e  >2 

2  2  J 

t-t  -iC.  (£)"  (58> 

2t  2  \  n  / 

where  £2-  {;  *  ic2^(S2/n)  is  small  and  imaginary. 

Equation  (57)  is  solved  numerically  and  in  Figure  5, 
the  variable  Re  [ (1-c+S)  (4/n82) **  is  plotted  as  a  function 

\c 

of  B/n  for  different  values  of  yfch  (y  *  0).  When 
y^-h  *  0,  the  curve  follows  Equation  (58)  for  very  small 
B/n*5,  and  decreases  to  zero  at  8/n*5  *  .118.  As  Y 
increases  form  zero,  the  curve  falls  more  rapidly  and  is  zero 
at  smaller  values  of  0/n  . 
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The  effect  on  linear  gain  of  finite  values  of  y 
is  shown  in  Figure  6. 


B.  Variable  Parameter  Wiggler 

For  a  variable  parameter  FEL  where  r  #  0,  it  is 
convenient  to  introduce  new  independent  and  dependent 
variables: 


q 


p+ Y 


th 


+  iY 


B 


A  e 


-ip2/2r 


Equation  (53)  can  then  be  rewritten  as  follows: 


3q 


B 


i  -iiaiilL 

=  -  —  (q-e  +  A  )  S(o,-u  )  e  2 

r* 

=  G(q) 


where 


2 (l-r+6 ) 

s2r 


7th  +  iY 


(59) 
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Figure  6.  Constant  Parameter  Wiggler  with  Frequency  Discrimination 
Re (6  )  vs  (62/n) 4  7 


If  the  exponential  is  approximated  by 


in/2rq‘  ,  _  m 


-  1  ss  —  — '  2  +  .  .  , 


2rq‘ 


Equation  (59)  simplifies  to: 


a2B 

3q2 


la 


+  A?B 


B 


*  G(q) 


(6 


where 


i  "  g2r2 


Numerical  integration  of  the  equation  without 
expansion  of  the  exponential  indicates  that  as  expected  for 
n/2q2r  <  1,  the  growth  rates  are  not  altered  substantially 
The  solution  of  Equation  (60)  is: 


B 


*q 


4i  A 


H(!| (A 
V 


•*“(1-1) 

,H f 


(A  q'  )q'*G(q'  ) 


pU-i) 

-  H0)  (A,q)  /  dq*H<2)  (Aq'j  q^Gtq') 


(6 


(1) 

where  H 

y 


(2) 

and  H 

y 


are  Hankel  functions  of  order 


y 


(4  +  iot, 


i 


As  p  tends  to  infinity  in  the  right  half  of  the  p-plane, 
1  2 

B-*-  —  exp  (^ip  /2V)  ,  and  consequently  A(p)  =  B  exp  (ip  /2T)  , 

P 

tends  to  zero. 

The  pulse  amplitude  is  obtained  from  the  inverse 
Laplace  transform 


a (v, -uq) 


■i  Tv  2/2  f  * 

~2^—J  dq  * 


(q-e)  v+  i 


Cq 


(q-£) 2 

2 

B  (q) 


1  >  v  >  0  .  (62) 


In  the  limit  of  an  infinitely  long  wiggler  where 
r  ■*  the  integral  in  Equation  (62)  may  be  evaluated  by 
deforming  the  contour  Cq  towards  the  Im  q  axis.  The 
integrand  is 
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2 


% 


4i  A 

1 


l 

.  (q  -  e)  .  / .  uir 

+ 1  2 —  iIm-t 


00  ( 1— i ) 

•J  dq'  Hy(1)(Aiq)q'  *G(q') 


q 


A i q  >>  1,  Imf  A i q)  >  0 


The  contour  approaches  a  saddle  point  at  q-e  ^  i  vr^  +  A  , 
and  the  contribution  of  this  saddle  point  increases 
exponentially  ^  expfA^F^v).  However,  if  a(v,-uQ),  (v<l) 
has  this  behavior,  it  will  not  be  possible  to  match  a  to 
the  exponentially  decaying  solution  ^  exp(-A  F  2  v),(v  >  1)  in 
order  to  satisfy  the  continuity  of  a  and  3a/3v  at  v  *  1. 
This  implies  that  the  boundary  condition  at  v  ■  1  can  be 
satisfied  in  the  limit  A  r!  +  «  only  if 

i 


“ (1-i) 


/ 


(») 


dq  ( A : q)  q2  G(q)  =  0 


(63) 


— 00  ( 1 — i ) 
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r- ^ 


where  the  contour  passes  above  the  origin.  Equation  (63)  is 
the  eigenvalue  equation  in  the  limit  A  T  -*■<*>. 

If  e  =  0,  the  integral  can  be  expressed  in  terms  of 
Whittaker  functions,  and  the  eigenvalue  equation  is  then 
given  by 


/«  +  <“]  e"E! 

—  OO  L  J 


=  e 


.  TT 

‘I 


-  1 
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Z 

o  3 


(Z  ) 

_  u  0 

4'2 


+  W 


i  (z  5 

1  jj  o 
4 ' 2 


(64) 


=  0 


where 
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-iA 


2 

1 


2 


and  the  contour  passes  above  the  origin. 
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In  the  limit  |Z  |  ^  |t|  >>  1, 

O  L 


W 


i  (Z  ) 
1  jJ  0 

4 ' 2 


\  (V2) 


,,  .*5  iiry/4  V, 

(2ir)  e  (2Zq)  6 


•i( 


.v, 


Ai  (2ZJ  (Z_+ip) 


where  K  is  the  modified  Bessel  function  of  order  p/2 

p/2 

and  Ai  is  the  Airy  function.  Thus,  the  solutions  of 
Equation  (64)  are  given  approximately  by  the  zeros  of  the 
Airy  function.  For  the  most  unstable  solution 


(Zo  +  ip) 
57 

(2Z  )  3 
0 


O  ** 

«■  *  J 


and  therefore 


7T 


<5„  a  2  e 
2 


(65) 
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Equation  (64)  is  solved  numerically,  and  in  Figure  7, 
the  variable  Re  ((l-r+<$)  (4/nB2)15)  is  plotted  as  a  function 
of  BT/n  .  As  sr/n  increases  from  zero,  the  curve 
decreases  from  (2)1*  and  asymptotes  to  zero  as  Br/n*5-*-00. 

The  effect  on  linear  gain  of  finite  values  of  yth 
and  y  has  been  calculated  by  numerical  solution  of 
Equation  (63) ,  and  the  results  are  displayed  in  Figures  7,  8 
and  9. 

A  more  exact  eigenvalue  equation  can  be  obtained  by 
replacing  the  Hankel  function  in  Equation  (63)  with 

the  appropriate  homogeneous  solution  of  Equation  (59),  namely 
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the  solution  B  +  which  asymptotes  to  %  exp(iA  q)  for 
h  1 

|q  |  >>  1 : 


»(l-i) 

J  dq  B*(q)  G (q)  =  0  .  (67) 

-oo  (1-i) 

Prom  numerical  solutions  of  this  equation,  the  linear  gain 

curves  relating  Re  I(l-r+6)  ( 4/n6 2 )  to  gr/n5*  is  plotted 

in  Figure  12  for  a  value  of  n/r  =1  (y  =  o  and 

th. 

y  3  0).  An  additional  parameter  n/r  is  necessary  to 
completely  describe  the  linear  gain  characteristics. 

However,  the  difference  in  the  eigenvalue  solutions  between 
the  more  exact  and  approximate  eigenvalue  equations  is  small 
when  n/r  <  1. 

C.  Combination  of  Constant  and 
Variable  Parameter  Wigqler 

In  the  constant  parameter  section  v^  >  v  >0,  the 

pulse  amplitude  a  is  given  by  Equation  (55)  which  may  be 

rewritten  in  the  form: 


a(v,-u  ) 
o 


a<o*-u0> 

j-lr^4 


(q ,-e) T^v 

6  ^  qj  (<*j“G+A:) 

“(q.-e) (2q.-e)-A*j 


(68) 
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where  q .  are  the  roots  of 


q2<q-e)2  -  Ajq2  +  i  a,  =  0 


In  the  variable  parameter  section  1  >  v  >  vi  , 
Equation  (47)  is  applicable  and  the  corresponding  equation 
for  B  is: 


32B 

w  +  A*B 


2 

e2r 


m 

2Tq2  .  _ 
e  -  1  }  B 


rU  1 


T"  Gi  +  (?)  e 


(69) 


where 


!i(q)  =  '  [^i7'*(v.’'u 


0)  +  (<3“E) 


G2(q) 


b  (vx)  end  b2  (Vj )  are  defined  by  Equations  (49a)  and 
(49b) . 


Thus,  the  eigenvalue  equation  analogous  to 
Equation  (61)  in  the  limit  Air*(l-v  )-*■  «  is 


dq  H  (Ajq)  q  e 


-» (1-i) 


(W  W 

v <*> +  2i°iKr^  +  -V" 


=  o  . 


When  e  «  0,  the  integral  can  be  expressed  in  terms  of 
Whittaker  functions. 

Equation  (70)  determines  as  a  function  of 

Oj,  r^v  ,  and  e* 

Equation  (70)  is  solved  numerically  and  in  Figure  10, 
the  curves  relating  Re  [  (l-r+6)  (4/nB2)**  ]  to  er/n*5  are 
plotted  for  different  values  of  rSr  . 


V.  SUMMARY  AND  DISCUSSION 


A  Summary  of  Linear  Gain  Results 

In  this  section,  the  linear  gain  characteristics 
per  pass  of  the  short  pulse  PEL  for  the  most  unstable  mode 
is  summarized.  However,  before  presenting  this  summary, 
the  relationship  between  the  symbols  introduced  in  the 
linear  eigenmode  analysis  and  the  physical  variables  of  the 
PEL  will  be  recapitulated. 


6  Gain  per  pass,  i.e.,  the 

amplitude  of  the  EM  pulse  on 
pass  n+1  is  e5  times  the  ampli¬ 
tude  after  the  nth  pass 


k 

w 


Higgler  wave  number 


k  A  =  B 
w  w  w 


L 


V 


k 

s 


Higgler  magnetic  field  amplitude 
Higgler  length 

Design  frequency  and  wave  number 
of  the  electromagnetic  pulse 


r  Fractional  reduction  per  round 

trip  of  pulse  amplitude  due  to 
losses  on  reflection  at  the 
mirrors 
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Resonant  energy  in  units  of  me2 


Total  number  of  electrons, 

/  ndl,  in  the  beam  per  unit 
area.  Note  that  the  beam  is 
considered  to  be  finite  and  uni¬ 
form  over  a  length  1  <<  L/2y2. 

at  is  the  pass-to-pass  temporal 
advance  of  the  electromagnetic 
pulse  relative  to  the  electron 
beam  on  entry  into  the  wiggler 


Ayr  is  the  change  in  y  of 
electrons  freely  accelerated  by 
the  accelerating  potential 


Aw  is  the  frequency  half  width 
of  the  frequency  discriminator, 
i.e.,  the  effective  reflection 

coefficient  is  r  (w) -t/  [1+i  (w-w)/Aw] 

s 


( Y )  -  is  the  difference 
between  the  mean  energy  of  the 
beam  and  the  resonant  energy 


Vth  is  the  mean  energy  spread 
of  the  electron  beam 


Fractional  length  of  the  initial 
constant  parameter  section  of  a 
variable  parameter  wiggler. 


.  No  Frequency  Discrimination 

For  an  FEL  without  frequency  discrimination,  the 
linear  gain  per  pass  Re  5  is  conveniently  expressed  in 
terms  of  the  dimensionless  variable 

i 

Re(i6^)  =  Re  (l  -  r  ♦  4  *  6Yth  )  (^§t) 

Figure  3  displays  a  plot  of  Re(i6^)  as  a  function 
of  the  parameter  6/n,  obtained  from  a  numerical  solution 
of  the  eigenvalue  Equation  (37)  for  a  constant  parameter 
wiggler.  Re(i£  )  is  zero  at  8/n  *  0.065  and  increases  in 
magnitude  as  6/n  decreases.  In  the  limit  of  a  very  long 
wiggler  where  6/n  -*■  0  while  62n  remains  finite,  the 
linear  gain  is 

Red  -  r  +  6  +  8yth)  -  |  (n&2)  1/3  •  (71) 

In  Figure  4,  Re(i6^)  is  plotted  as  function  of 
28r  Vn  for  a  variable  parameter  wiggler.  The  graph  is 

obtained  from  a  numerical  solution  of  Equation  (43). 

V2 

Re(i6.)  is  zero  at  28T  /n  ~  1.42,  and  increases  in 
1  % 

magnitude  as  28T  /n  decreases.  In  the  limit  of 

3/2 

28T  /n  -*•  0,  the  linear  gain  is 

3  3/2  V3 

R(1  -  r  +  5  +  B7th>  j  (nB2)  .  (72) 
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The  linear  gain  of  a  variable  parameter  wiggler  can 
be  increased  by  the  addition  of  a  constant  parameter 
wiggler  section  at  the  front.  In  Figure  4,  the 

3/, 

corresponding  curve  of  Re(i  5^  )  as  a  function  of  26T  /n 

is  plotted  for  a  composite  constant  and  variable  parameter 

wiggler.  It  is  obtained  from  a  numerical  solution  of 

Equation  (50),  assuming  that  the  fractional  length  v  <<  1 

of  the  constant  parameter  section  is  given  by  F^v1=  1.5. 

The  enhancement  of  linear  gain  increases  with  F^v  , 

3/  1 

and  can  be  particularly  significant  for  2gr  2/n  2  1.42. 

Deviations  y  of  the  mean  energy  of  the  beam  from 
the  resonant  energy  do  not  affect  linear  gain.  However,  a 
finite  energy  spread  y  decreases  linear  gain  by  gy^. 

2.  Frequency  Discrimination 

For  a  FEL  with  frequency  discrimination  (v  =  1/8) 

the  linear  gain  per  pass  Re  5  is  conveniently  expressed  in 

_ 

terms  of  the  dimensionless  variable  Re<$2  *  Re[  (l-r+6)  (4/nB  )  ]. 

In  Figure  5,  the  constant  parameter  linear  gain 
curve  (obtained  from  a  numerical  solution  of  Equation 
(57)]  which  relates  Re  S2  to  the  parameter  6/n 

is  plotted  for  different  values  of  the  energy  spread  Yth 

h 

(but  with  y  *  0) .  Linear  gain  is  possible  when  8/n 

is  below  an  upper  limit  (8/n15)  .  For  a  beam  with  7=0, 

max 

(S/n*5)  «  0.118.  In  the  limit  8/n 2  -  o,  the  linear 

max 

gain  is: 
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Re  (1 


r  +  5) 


(73) 


_  ^ 

As  Yth  increases,  (B/n  )max  decreases,  and  the 
magnitude  of  Re  <$2  is  reduced. 

In  Figure  7,  the  variable  parameter  linear  gain 
curve  [Equation  (63)]  which  relates  Re  62  to  the 
parameter  Br/n2  is  plotted  for  different  values  of 

—  k  _  — 

Yth/r  (but  with  Y  ■  0).  For  a  beam  with  Yfch  *  0 

Re  ( 1  -  r  +  6 )  -*• 

A  finite  energy  spread  y  .  reduces  the  magnitude  of  Re  <5  . 

tn  z 

In  Figure  10,  the  linear  gain  curve  of  a  composite 

constant  and  variable  parameter  wiggler  [Equation  (70)] 

which  relates  Re  52  to  Bryn*5  is  plotte  I  for  different 

values  of  v^.  The  addition  of  a  constant  parameter 

section  can  enhance  the  linear  gain  of  a  variable  parameter 

H 

wiggler  and  can  be  of  significance  for  BT/n  2  1.0. 

For  a  given  set  of  FEL  parameters,  the  linear  gain 
is  not  in  general  a  maximum  at  y  *  0.  The  effect  on  the 
linear  gain  of  a  variation  in  y  is  shown  in  Figure  6  for 


(nS 2/2) 


BT 

* 


V  +  0 


BT 

r\h 


(74) 
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a  constant  parameter  wiggler,  and  in  Figure  9  for  variable 
parameter  wiggler s.  Some  enhancement  of  linear  gain  is 
possible  by  making  an  appropriate  non-zero  choice  of  y  . 

These  results  are  applicable  to  the  special  case  of 
v  3  1/6.  The  rationale  for  studying  this  case  is  based  on 
the  observation  that  frequency  discrimination  has  an  effect 
(proportional  to  1/v)  which  tends  to  retard  the  pulse 
advance,  that  is,  it  tends  to  reduce  the  effective  pulse 
advance  produced  by  a  finite  positive  value  of  6.  To 
avoid  the  EM  pulse  running  either  ahead  of  or  lagging 
behind  the  electron  pulse  it  is  indicated  that  one  should 
choose  v  =  1/6.  Some  numerical  evidence  from  our 
simulations  indicates  that  this  choice  is  indeed  nearly 
optimal. 

The  graphs  which  have  so  far  been  described  are 
obtained  from  approximate  eigenvalue  equations  derived 
after  expansion  of  the  exponential  term 

(• 


31 

■m  p 

2  3p 


-  1/  =  -  ^ 


31 

in  p 

2  3p 


in  Equation  (30)  and  Equation  (53)  . 

This  expansion  of  the  exponential  should  be 
reasonable  for  the  constant  parameter  wiggler  when 
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V  w 

(nB2)  3  <<  1  and  ( n 6 2 >  <<  1,  without  and  with 

frequency  discrimination  respectively.  The  dominant  poles 

of  the  integrand  in  the  integral  of  Equation  (34)  [Equation 

(54)]  are  then  determined  to  a  good  approximation  by  the 

roots  of  the  cubic  Equation  (38)  [quartic  Equation  (56)]. 

In  order  to  assess  the  error  introduced  for 

the  variable  parameter  wiggler  without  frequency 

discrimination,  the  more  exact  eigenvalue  Equation  (42)  is 

solved  numerically,  and  the  linear  gain  curve  relating  Re(i6 

3/  1 

to  the  parameter  26T  2/n  is  plotted  in  Figure  11  for  a 

value  of  n/r  *>  1.0  (Y  *  0,  Y..  =  0).  One  additional 

parameter,  namely  n/r,  is  necessary  to  completely 

describe  the  linear  gain  characteristics.  Similarly,  from 

the  more  exact  eigenvalue  Equation  (67)  for  the  variable 

parameter  wiggler  with  frequency  discrimination,  the  linear 

\c 

gain  curve  relating  Re  «2  to  the  parameter  BT/n  is 
plotted  in  Figure  12  for  a  value  of  n/r  *  l.o  (y=0,  Yfch=  0)  . 

In  both  the  FEL  with  and  without  frequency 
discrimination,  the  difference  in  eigenvalue  solutions 
between  the  more  exact  and  approximate  eigenvalue  equation 
is  small  when  n/r  <  1. 

2.  Significance  of  Results  for  FEL  Design 

The  parameter  regime  of  interest  for  variable 
parameter  FEL's  is  dictated  not  only  by  the  need  to  grow 
the  electromagnetic  pulse  to  a  large  amplitude  in  a  finite 
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Figure  11.  Variable  Parameter  Wiggler  without  Frequency 

3/ 

Discrimination  -  Re(i6  )  vs  2gr  2/n 


number  of  passes  (i.e.,  linear  growth  characteristics),  but 
also  by  the  requirements  of  effective  electron  trapping  and 
stable  propagation  of  the  finite  amplitude  pulse.  Although 
the  latter  two  requirements  are  not  the  main  topics  of  this 
paper,  it  will  be  useful  to  briefly  indicate  some  con¬ 
straints  on  the  choice  of  parameters  which  they  impose. 

The  comments  on  the  constraints  imposed  by  the  requirements 
of  non-linear  propagation  will  be  abstracted  from  Reference 
1.  They  are  compared  to  those  summarized  for  linear  gain. 
The  electron  equations  of  motion  [Equations  (14)  and  (15)] 
are  derivable  from  the  Hamiltonian  H  : 


H 


li  . 


I>  -  a  cos  (ip  +  £) 
s 


Y  +  v(^) 


In  a  finite  amplitude  stationary  state  where  a  and  C 

s 

are  constant,  electrons  with  small  Y  are  trapped  in  the 
potential  well  v(<ji).  The  effective  potential  well  depth 
is 


AV  =  2a  cos ^ 

s  r  r 


TT . 

smij>r 


where  ^  <  j 


is  given  by 


a  sin  ip  =  r 
s  r 

The  oscillation  'frequency'  at  the  bottom  of  the  potential 
well  is 


(a  cos  4*  ) 
s  r 

The  magnitude  o£  the  stationary  value  of  a  may 
be  estimated  from  energy  balance.  Assuming  ag 
independent  of  v,  Equation  (52)  yields 


(1-r)  ag2  x  n  ag  sin  (4;  +  £)^> 


n 


r 


where  ffc  is  the  fraction  of  trapped  electrons,  which 
should  be  close  to  unity. 

From  these  properties  of  the  potential  well,  the 
following  observations  can  be  made: 

1)  To  secure  a  finite  well  depth,  it  is  desirable 
to  have: 


sin  4/ 

r 


A 

a 

s 


(75) 
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2)  To  trap  a  large  fraction  of  electrons,  the 


energy  spread  y  should  satisfy 


< 


(76) 


We  will  discuss  below  how  this  nonlinear  limitation  on 
thermal  spread  compares  with  that  determined  by  the 
limiting  spread  for  linear  gain  calculated  herein. 

In  reference  1,  it  is  shown  that  the  gain  of 
sideband  instabilities  driven  by  electrons  trapped  in  the 
ponderomotive  potential  wells  which  oscillate  at  frequency 
(ag)?s  is  of  the  order  of  2-3  times  the  gain  experienced 
by  the  signal  (which  is  by  definition  equal  to  the  losses 
in  the  reflectors  in  steady  state) .  Furthermore,  sideband 
frequencies  in  the  range  0.2(ae)*  to  (ae)^  show  such  gain. 

s  s 

We  can  therefore  make  a  crude  estimate  of  the  needed 
frequency  discrimination. 

To  suppress  sideband  frequencies  of  the  order  of  a 

1  W 

quarter  of  the  oscillation  frequency  -  -•  (a  )  *,  the 
reduction  in  power  due  to  filtering  by  frequency  discrimi¬ 
nation  should  exceed  the  relative  gain  due  to  sideband 
instabilities: 
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r 


_i_ 

16 


a  _ 

>  2  (1  -  r)  (G  -  1) 


G  is  the  ratio  of  gain  of  the  sideband  instability  to  the 
gain  of  the  pulse  amplitude.  Substituting  for  a  ,  and 

S 

taking  v  *  -g-s 

g2r2  >  32 (G-l) 

n  fth 


(l-r) r 
n 


Such  arguments  are,  of  course,  not  quantitative, 
but  the  results  of  our  numerical  simulation  suggest  that  a 
rough  rule  compatible  with  the  above  discussion  is 


B2 


rn 

(l-r)  3 


*5 


70  , 


(77) 


where  the  numerical  coefficient  is  uncertain  by  a  factor  of 
the  order  of  2. 

We  note,  of  course,  that  large  6  (i.e.,  high 
frequency  discrimation)  is  desirable  for  stable  propagation 
in  the  presence  of  the  sideband  instability.  Comparing  now 
with  the  linear  theory  of  this  paper,  we  note  that  Figure  6 
gives  a  plot  of  (l-r+5)  (4/n@z)is  vs.  QT/rfi.  Here,  of 
course,  small  8  is  desirable  for  linear  gain.  At 
marginal  stability  (6*0)  and  using  the  lowest  value  of  B 
allowable  from  Equation  (77),  i.e.,  6*8.5 (l-r)3/rn)  *  we  see 
that  Figure  6  may  be  regarded  as  a  transcendental  equation 
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determining  the  maximum  value  of  (l-r)r/n  for  which  a 
value  of  frequency  discrimination  may  be  chosen  which  both 
allows  for  linear  gain  and  stable  propagation.  Solving 
this  equation  numerically/  we  find 

(1  -  r)  r/n  <  .05  (78) 

as  a  fundamental  limitation  on  the  design  of  a  variable 
parameter  short  pulse  FGL. 

Obviously  to  provide  for  reasonable  positive  gain, 
we  must  make  (l-r)r/n  somewhat  less  than  this  critical 
value.  The  situation  may  be  ameliorated  somewhat  by  use  of 
a  constant  wiggler  section  (see  Figure  6)  to  allow  linear 
gain  at  higher  8.  However,  the  use  of  such  a  section 
leads  to  saturation  of  the  linear  growth  phase  at  an 
amplitude  which  may  be  too  low  for  trapping  since 

a  m 

( a sat>  V!  «  it.  It  appears  then  that  in  practice  the 
length  of  the  constant  wiggler  should  be  limited  by 
r*v  <  i.  with  this  optimization,  Equation  (78)  appears 
as  a  reasonable  design  limitation.  Note  that  from 
Equation  (75)  this  limitation  may  also  be  expressed  in  the 
form 


sin  <Pr  <  0.25 
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We  note  also  that  using  this  value,  Equation  (76) 
foe  the  nonlinear  limitation  on  temperature  may  be 
expressed  as: 


W  <  2 


while  from  Figure  (8)  we  see  that  the  linear  gain  is 

_  k 

strongly  reduced  if  y  h/r  >  1.  Hence,  we  see  that  good 
linear  gain  requires  about  twice  as  small  an  effective 
energy  spread  as  would  be  deduced  from  looking  at  nonlinear 
trapping  alone. 

We  have  seen  that  the  properties  of  a  very  short 
pulse  FEL  can  be  characterized  by  the  dimensionless 
parameters  r  (related  to  the  total  variation  of  resonant 
energy  induced  by  the  variable  parameters) ,  n  (related  to 
the  strength  of  electron  current)  and  &  the  advance  of 
the  electromagnetic  pulse  compared  to  the  next  pulse  in  the 
electron  train.  We  have  also  introduced  a  frequency  dis¬ 
crimination  v  in  order  to  suppress  sideband  instabilities 
in  the  nonlinear  evolution.  In  terms  of  these  parameters, 
we  have  calculated  the  linear  gain  of  the  system  and  its 
degradation  by  thermal  spread.  We  note  that  small  thermal 
spread,  low  frequency  discrimination,  high  current,  and 
small  variation  of  resonant  energy  are  conducive  to  good 
linear  gain. 
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In  our  companion  paper,  we  will  discover  that  large 
frequency  discrimination  is  required  to  suppress  sideband 
instabilities  in  the  nonlinear  region,  thus  determining  a 
minimum  value  of  n/T  for  which  good  operation  may  be 
expected  with  optimal  frequency  discrimination  if  beam 
energy  spreads  and  emittance  can  be  kept  low  enough. 
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Simulation  of  Short  Electron  Pulse  Free  Electron 
Lasers  with  Variable  Parameter  Higglers 

B.  N.  Moore,  M.  N.  Rosenbluth, 
and  H.  Vernon  Wong 

I .  INTRODUCTION 

The  free  electron  laser  (FEL)  is  a  device  for 

converting  the  directed  energy  of  a  relativistic  electron 

beam  into  short  wavelength  electromagnetic  energy.  This  is 

achieved  by  exploiting  a  resonant  interaction  between  the 

beam  and  wave  in  a  region  where  beam  dynamics  is  dominated 

by  a  spatially  varying  magnetic  field.  '  Such  devices 

may  be  operated  in  either  an  simplifier  or  oscillator  mode 

and  experimental  lasers  have  been  operated  in  each  mode.^ 

Analyses  of  constant  parameter  wiggler  FELs  have  identified 

2 

the  major  limitations  on  energy  extraction  efficiency. 

The  most  important  is  that  the  EM  pulse  tends  to  grow  to 
a  level  at  which  electrons  just  become  trapped,  but  because 
additional  extraction  would  bring  electrons  below  the 
resonant  energy,  further  growth  is  prevented.  The 
saturated  state  then  is  one  in  which  efficient  energy 
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extraction  takes  place  only  near  the  end  of  the  wiggler. 

The  variable  parameter  wiggler  is  a  modification  of  the 
2 

FEL  which  overcomes  this  problem  by  optimally  adjusting 
the  resonant  energy  as  a  function  of  position  through  the 
introduction  of  suitable  wiggler  amplitude  and/or  wave¬ 
length  variation.  Alternatively,  one  could  in  principle 
employ  an  electrostatic  field  to  make  up  for  radiative 
losses  and  keep  an  electron  in  resonance  with  a  constant 
wiggler.  It  was  anticipated  that  the  variable  parameter 
wiggler  would  have  problems  associated  with  trapped  parti¬ 
cle  instabilities  at  sidebands  offset  from  the  resonant 
frequency.  However,  these  instabilities  can  be  suppressed 
by  the  introduction  of  frequency  discrimination  into  the 
FEL  optics.  Such  frequency  discrimination,  of  course, 
reduces  growth  in  the  unsaturated  state.  One  of  the 
accomplishments  of  this  simulation  effort  has  been  to 
demonstrate  that  efficient  energy  extraction  may  be 
expected  in  interesting  parameter  regimes  with  frequency 
discrimination. 

The  focus  of  this  simulation  study  has  been  the  FEL 
with  variable  parameter  wiggler.  A  family  of  simulation 
codes  has  been  written  and  applied  to  an  investigation  of 
growth  and  saturation  of  FELs  with  variable  parameter 
wigglers.  This  report  describes  the  codes,  the  results 
obtained  to  date  and  future  plans.  Some  correlation  of 
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simulation  results  with  other  analyses  and  experiments  has 

3  4 

been  presented  elsewhere.  '  This  report  is  organized 
into  sections  describing  the  physical  and  mathematical 
model  of  an  FEL  with  discussion  of  the  underlying  assump¬ 
tions  and  simplifications.  This  is  followed  by  descrip¬ 
tions  of  the  numerical  algorithms  including  considerations 
of  numerical  stability  and  error  estimation. 

Due  to  computer  time  limitations,  most  of  the 
simulations  have  been  done  using  ultrashort  electron  pulses 
modeled  by  a  6-function.  The  results  should  apply  to  cases 
where  the  electron  pulse  is  much  shorter  than  the  slippage 
distance;  i.e.,  the  difference  in  distance  traveled  by 
photons  and  electrons  in  the  time  taken  by  a  photon  to 
traverse  the  wiggler. 

The  results  of  the  simulation  of  various  cases  are 
summarized  in  the  following  sections  and  particularly 
interesting  simulations  are  presented  in  greater  detail. 

In  a  separate  report"*  we  have  made  extensive  analytic 
studies  of  linear  growth  of  FEL  oscillators,  and  we  will 
present  in  this  report  some  confirmations  of  those  studies. 
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II.  FEL  PHYSICS 


The  essential  features  of  an  FEL  oscillator  are 
shown  schematically  in  Figure  1. 

It  consists  of  an  EM  cavity,  bounded  on  each 
end  by  mirrors,  into  which  are  injected  pulses  of  high 
energy  electrons.  The  mirrors  also  serve  to  extract 
EM  energy  and  possibly  as  filters.  Also  within  the 
cavity  is  a  region  of  spatially  varying  magnetic  field  with 
its  major  component  transverse  to  the  direction  of  pulse 
propagation.  This  is  the  wiggler,  whose  amplitude  and 
wavelength  may  be  adjusted  to  optimize  energy  extraction. 
The  electron  pulses  and  coherent  EM  pulses  pass  synchro¬ 
nously  through  the  wiggler  field.  The  effect  of  the 
wiggler  field  acting  alone  would  be  to  establish  the 
transverse  motion  of  the  electrons.  In  fact  the  EM  field 
couples  with  the  magnetic  field  of  the  wiggler  to 
establish  an  effective  ponderomotive  potential  for  the 
electrons.  Both  fields  are  assumed  circularly  polarized. 

The  electron  equation  of  motion  in  terms  of  energy 
Y  and  phase  angle  are  given  below.  A  full  derivation 

may  be  found  in  Reference  2. 
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Where 


ip  is  particle  phase 

?  is  electromagnetic  phase 

A  is  electromagnetic  vector  potential  amplitude 

s 

A  is  wiggler  vector  potential  amplitude 
w 

<|)  is  the  electrostatic  potential  which  mocks 
0  up  the  wiggler  parameter  variation. 

X  is  the  EM  wave  length 

9 

Xw  is  the  wiggler  wave  length 
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m,  e  and  c  are  electron  mass,  electron  charge  and  light 
speed,  y  is  the  relativistic  factor  for  the  beam. 


Y 


Electron  Energy 
me2 


(6) 


The  transverse  electron  current  serves  as  a  source 
for  the  electromagnetic  field  of  the  EM  pulse  governed  by 
Maxwell 1 s  equations . 
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The  current  densities  may  be  calculated  from  the  electron 
trajectories 
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Where  the  notation  ^  ^  implies  averages  over  the  initial 
phase  and  energy  distribution  of  the  electron  beam,  and  the 
function  n  is  the  number  density  of  the  electron  beamlet, 
and  vr  is  the  resonant  y.  Here  it  has  been  assumed  that 
the  EM  field  is  of  the  form 


2  1  3A 

E  -  -  ^ 

c  3t 


(11) 


eA 


me 


a  (z,t)  x  cos  (ip  )  -  y  sin  (ip  ) 
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(14) 


w  =  k  c 
s  s 


(15) 


where  and  ?  are  slowly  varying  functions  of  z  and  t 
and  are  independent  of  x  and  y. 

The  EM  pulse  leaves  the  wiggler,  propagates  to  the 
mirror  where  it  is  filtered  and  partially  reflected  before 
encountering  the  next  electron  pulse.  These  processes  are 
modeled  as  follows. 


a 

3  into  filter 
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(reflection) 
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Where 

o>  is  the  bandpass  cutoff  frequency  centered 
F 

about  (j 

s 

r  is  the  extraction  mirror  reflection  coefficient 

I*  is  the  wiggler  length. 

Before  entrance  into  the  wiggler  on  the  next  pass, 
a  fixed  time  delay  is  introduced.  The  electron  beam 
density  into  the  wiggler  is  assumed  to  be  perfectly 
periodic,  with  the  period  At. 

The  major  assumptions  underlying  this  model  are: 

(1)  The  wiggler  field  has  negligible  axial  and 
fringe  components. 

(2)  Transverse  variations  of  all  quantities  are 
neglected . 

(3)  An  EM  pulse  is  a  nearly  plane  wave  whose 
amplitude  and  phase  vary  negligibly  over  a  few 
EM  wavelengths. 

(4)  The  appropriate  EM  filter  with  Lorentzian  line 
shape  is  assumed. 

(5)  y  >>  1. 
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III.  SIMULATION  EQUATIONS 


By  appropriate  normalization,  the  FEL  equations  may 
be  further  simplified  for  the  simulation  code.  Other 
simplifying  assumptions  will  also  be  convenient. 
Additionally,  finite  difference  algorithms  must  be  selected 
which  do  not  introduce  numerical  instabilities  and  which 
efficiently  approximate  the  time  and  spatial  derivatives 
involved. 

Introduce  the  variables  u  and  v. 


u 


(18) 


(19) 


u  ■  constant  follows  an  electron,  and  v  *  constant 
a  photon.  Lw  is  wiggler  length,  V  the  electron 
longitudinal  velocity.  The  electron  beam  runs  from  -uQ  to 
zero.  The  wiggler  interaction  occurs  for  1  >  u  +  v  >0 
(See  Figure  2).  The  equations  for  this  region  of  inter¬ 
action  are: 


2 


Figure  2.  Simulation  Coordinate  System 
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where  a,  r#  n  and  h(u),  expressed  in  terms  of  physical 
parameters,  are  defined  by 
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Z  is  nominal  electron  pulse  length 

r  is  electron  beam  radius 

o 

I  is  the  beam  current 


In  terms  of  the  normalized  variables,  the  pass-to-pass 
EM  processes  are  modeled  as: 


3a<V)  (n\ 

-  +  v  a_  (v)  =  v  a  (u  =  0 ,  v) 

9v  F 


.  (n+1) 
a 


(u 


u  ,  v)  =  ra_ (v+s) 

0  r 


k  L 
w  w 

k  c 
s 


U> 


F 


r  =  reflection  coefficient  of  extractor 


k  c 

fi  =  - — - ?  At  .  =  : — — - —  At 


Where 


with  At  being  the  time  delay  of  the  electron  pulse 
rel 

with  respect  to  the  EM  pulse  at  the  wiggler. 

Some  preliminary  attempts  at  algorithms  for  solving 
these  equations  met  with  limited  success,  either  due  to 
problems  with  numerical  instabilities  or  lack  of 
computational  efficiency.  The  EM  amplitude  and  phase 
(a  .5)  and  the  particle  phase  space  variables  (y,<M  are 
defined  on  a  discrete  grid  in  the  u  -  v  plane.  The 
indices  i  and  i  define  the  location  on  this  grid. 

The  index  i  is  the  particle  label.  The  following 
algorithms  were  eventually  selected. 

u  =  u  +  (i  -1)  Au  (36 ) 

o  u 

v  *  v  +  (i  -  1)  Av  (37 ) 

0  V 

Y  .  (u,v  +  Av)  =  y  .  (u,v)  -  Av  (a  sin  'k  +  a  cos  ^  .  J 
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(38) 

Where  the  centered  EM  field  components  ares 
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The  particle  phases  are  advanced  for  a  fixed  u  value. 
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The  contribution  of  the  wiggler  parameter  variation  to  the 
particle  phases  is  isolated  to  the  term: 
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Steps  (38)  through  (45)  are  repeated  for  each  particle. 

Then  the  particle  currents  are  calculated  and  the 
field  amplitudes  at  the  next  vaue  of  u  are  determined. 
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Steps  (46)  through  (51)  are  repeated  for  each  v  interval 
in  the  wiggler  (0  <  u  +  u  <  1).  For  other  values  of  v 
these  quantities  propagate  unchanged.  Finally,  the  entire 
series,  (36)  through  (51),  is  repeated  for  each  value  of  u. 
The  basic  scheme  is  depicted  in  Figure  2.  Specification  of 
the  EM  pulse  as  a  function  of  time  into  the  wiggler  is 
equivalent  to  specification  vs  v  at  u  «  “u0* 

Similarly,  the  electron  pulse  density  distribution 
is  a  specified  function  h(u)  at  v  *  -  u.  The  field 
equations  are  integrated  along  the  vertical  characteristics 
(vaconstant)  and  the  equations  of  motion  along  horizontal 
(u*constant)  characteristics.  Mote  that  the  y,  sin  ij/  and 
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cos  functions  are  advanced  independently  with  an 
algorithm  which  preserves  sin2  4>  +  cos2  ij.1  .  The  particles 


at  a  given  u  are  advanced  from  v  *  -u  to  v  ■  1-u  with 
an  algorithm  which  estimates  the  effects  of  the  averaged 
EM  field  acting  on  the  particles.  Once  the  phase-averaged 
particle  currents  are  available  at  a  given  u,  the  EM 
fields  are  advanced  to  the  next  u  value.  The  process  is 
repeated  until  the  electron  pulse  and  EM  pulse  have  moved 
completely  through  the  interaction  region.  The  EM  pulse 
then  propagates  undisturbed  to  the  mirrors  where  energy 
extraction  and  optical  discrimination  are  carried  out. 
Negligible  interaction  on  the  return  through  the  wiggler  is 
assumed. 

One  group  of  simulation  programs  incorporates  the 
additional  renormalizations 
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This  renormalization  involves  a  reduction  of  the 
output  amplitude  a  to  a  such  that  the  input  pulse 
energy  is  always  at  a  level  aQ2 ,  small  enough  to  be  in  the 
linear  regime. 


This  keeps  the  EM  pulse  at  a  constant  total  energy 
before  entry  into  the  wiggler.  The  result  of  many  repeti¬ 
tions  of  this  algorithm  will  be  a  pulse  corresponding  to 
the  fastest  growing  linear  eigenmode,  which  can  be  compared 
with  analytic  results  and  used  to  begin  the  nonlinear 
simulation. 

Three  options  are  provided  for  loading  the  initial 
electron  beam.  These  are:  cold  beam  and  Maxwellian  and 
Lorentzian  thermal  distributions.  The  distributions  are 
specified  by  ^  y  ^  ,  which  is  the  mean  y,  and 
which  characterizes  the  thermal  spread. 
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The  particles  are  also  distributed  uniformly  in 


phase . 
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IV.  DIAGNOSTICS 


A  number  of  diagnostics  are  provided  in  the  codes, 
including  growth  rates,  integrated  power,  phase  space  plots 
and  Fourier  transforms.  The  growth  rate  represents  the 
pass  to  pass  exponentiation  as  follows: 


Re  6 


(60) 


The  particle  phase  space  corresponding  to  any  point 
in  the  interaction  region  me\y  be  displayed  for  any  pass, 
allowing  a  determination  of  trapping  efficiency.  The 
energy  extracted  per  particle  is  monitored  as  a  function 
of  pass  with  the  parameter: 


N, 


5y 


ext 


E  E  (v-/ r<sv)  h(u)  ^ 

in  io=l  '  v,  ' 


(61) 


where  Ne  is  the  total  number  of  simulation  electrons. 

The  Fourier  transform  diagnostics  and  their 
relation  to  physical  parameters  are: 
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a  +  i  a 
c  s 


(62) 


a  (v)  = 

Ak  =  £  *(v) 
iv 


i< (ivAv) 
e  v 


K 


k  L 


w  w 


(V/c) 


(63) 

(64) 


Note  that  the  diagnostic  frequency  (or  wave  vector)  is 
expressed  relative  to  the  nominal  resonant  wavelength. 
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V .  RESULTS 


The  major  tasks  accomplished  to  date  include  the 
verification  of  the  code  against  established  theory  of  the 
constant  wiggler  PEL,  confirmation  and  extension  of  the 
linear  theory  of  variable  wiggler  devices,  including 
thermal  effects,  and  investigations  of  the  nonlinear 
phenomena  associated  with  saturation  of  both  short  and  long 
pulse  FELs.  These  investigations  have  established  the 
necessity  for  some  frequency  discrimination  in  order  to 
achieve  acceptible  performance  for  both  classes  of  devices, 
and  the  scaling  laws  for  variable  parameter  wiggler  FELs 
with  frequency  discrimination  have  been  determined. 

Preliminary  simulations  of  finite  length  electron 
beams  without  frequency  discrimination  have  given  results 
similar  to  those  of  Goldstein,  et  al.^  for  parameters 
associated  with  the  LANL  proposed  experiments.  In  the 
absence  of  frequency  discrimination,  the  EM  pulse  shape  at 
saturation  is  irregular  and  characterized  by  random  peaks. 
With  frequency  discrimination,  much  smoother  pulse  shapes 
were  obtained,  and  the  saturated  pulse  energy  was  somewhat 
smaller.  However,  the  variation  with  pulse  length 
indicates  that  with  longer  beam  pulses  such  as  those 
actually  envisaged,  frequency  discrimination  would  be  more 
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effective  in  producing  not  only  smoother  saturated  pulse 
shapes,  but  also  higher  pulse  energies.  Further  work  on 
finite  length  pulse  simulation  will  be  published 
separately. 

In  the  following  sections,  we  focus  exclusively  on 
the  limit  of  very  short  electron  beams  in  which  the 
electron  pulse  length  l  is  short  compared  to  the  slippage 
distance  ^Ly/kg  >:>  ^  Jl  or  r^|uQ|  <<  1. 

In  the  simulation,  this  limit  is  modeled 
as  a  6-function;  that  is,  h(u)  ■  5(u).  Obviously 
this  model  requires  minimal  simulation  time  (only  one 
value  of  u)  so  that  long,  high  extraction  wigglers  may  be 
studied. 

As  a  check  on  the  accuracy  of  the  code,  the  linear 
gain  per  pass  calculated  from  the  linear  eigenmode  analysis 
was  compared  with  that  observed  in  the  simulations  for  a 
range  of  FEL  parameters.  In  the  parameter  regions  where 
the  theory  is  applicable,  the  agreement  between  theory  and 
simulation  is  satisfactory  (see  Figures  3,  4,  8,  9  and  10) 
and  is  of  the  order  of  or  better  than  15%. 

1.  No  Frequency  Discrimination 

In  FELs  without  frequency  discrimination,  the 
linear  gain  per  pass  for  a  very  long  constant  parameter 
wiggler  ( 8/ n  •+  0)  is  given  by 


Re  <5 

1 


Re 


r  i 


1  -  r  +  6  +  S  Yth 


(65) 


If  the  beam  and  wiggler  parameters  are  kept 
constant.  Re  |  (6/nj^3  j  increases  from  zero  as  the  phase 
advance  8  of  the  EM  pulse  is  increased,  reaches  a  maximum 
Re  |  (  8/n  ) <5  ij  ■  .047  at  B/n  *  .022,  and  finally 
decreases  to  zero  at  8/n  ■  .065.  The  theoretical  linear 
gain  curve ^  is  plotted  as  a  solid  line  in  Figure  3, 
while  the  linear  gain  observed  in  the  simulations  is 
represented  by  the  circled  points. 

If  the  wiggler  is  made  variable  by  applying  a 
finite  accelerating  electric  field  (T  4  0),  the  linear  gain 

remains  close  to  that  of  the  constant  parameter  wiggler 

3/ 

when  26  r  2/n  <<  1*  (The  accelerating  electric  field 
produces  essentially  the  same  effect  as  varying  the  wiggler 
parameters,  i.e.,  keeping  trapped  electrons  in  resonance  as 
they  radiate  energy) .  As  the  magnitude  of  the  variation  is 
further  increased,  Re  <5  j  decreases.  In  the  limit  of  a  very 
long  variable  wiggler  where  B/n  <<  i  and  r  >>  1, 

Re  «  0  when  2BT^/n  *  1.42.  In  Figure  4,  the 
theoretical  linear  gain  curve  is  plotted  as  a  solid 
line  indicating  the  variation  of  Re  6]  as  a  function  of 
26  r^/n.  The  results  of  the  simulations  are  represented  as 
circled  points  for  values  of  6/ n  “  3.3  x  10"^,  .01. 
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Without 


The  differences  between  the  observed  linear  gain  and  the 
theoretical  linear  gain  curve  for  S/n  *  *01  are  due  to 
finite  length  effects  not  included  in  the  theory.  For  most 
cases  of  physical  interest,  the  criterion  6/S  >  1  for 
neglecting  finite  length  effects  is  well  satisfied. 

The  pulse  shape  of  the  linear  eigenmodes  are 
displayed  in  Figure  5.  For  a  constant  parameter  wiggler, 
the  pulse  amplitude  increases  smoothly  to  a  maximum  and 
then  falls  to  zero  at  the  back.  As  the  wiggler  is  made 
variable,  the  change  in  the  pulse  phase  angle  in  going  from 
the  front  of  the  pulse  to  the  back  progressively  increases. 
At  the  same  time,  the  pulse  amplitude  becomes  more  sharply 
peaked  at  the  back. 

In  a  typical  run  in  which  the  EM  pulse  grows  from 
noise  levels  to  saturation,  the  early  phase  of  the 
evaluation  is  characterized  by  growth  of  the  linear 
eigenmode.  Eventually,  nonlinearities  reduce  and  then 
saturate  wave  growth. 

For  constant  parameter  wigglers,  the  maximum  pulse 
energy  occurs  for  values  of  the  pulse  phase  advance  parame¬ 
ter  0  which  are  small  and  positive.  This  decreases  as  S 
increases.  Similar  results  were  previously  obtained  by 
Colson  for  finite  length  electron  beams. ^  The  effect  of 
B  >  0  is  to  shift  the  EM  pulse  ahead  of  the  electron  beam 
on  each  pass  through  the  wiggler.  As  0  increases,  the 
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16 


Figure  5b.  EH  Pulse  Phase  vs.  v 


Figure  5d.  EM  Pulse  Phase  vs.  v 


net  interaction  time  of  each  section  of  the  pulse  with  the 
electron  beam  decreases  and  this  leads  to  a  progressive 
reduction  of  the  saturated  pulse  energy. 

For  a  variable  parameter  wiggler,  the  evolution 
enters  a  nonlinear  phase  during  which  the  peak  in  the  pulse 
amplitude  moves  to  the  front,  which  continues  to  grow  at 
approximately  the  linear  rate  after  the  back  has 
saturated.  Growth  continues  until  the  peak  amplitude  at 
the  front  of  the  pulse  is  large  enough  to  begin  to  trap  the 
beam  electrons.  At  saturation,  without  optical  discrimina¬ 
tion,  the  pulse  remains  peaked  at  the  front  and  never 
broadens  towards  the  back.  Most  of  the  energy  resides  in 
the  front  part  of  the  pulse  which  no  longer  interacts 
coherently  with  the  beam. 

The  pulse  shape  and  pulse  spectrum  at  saturation 
for  a  representative  simulation  run  with  the  following  FEL 
parameters  are  shown  in  Figure  6: 

-  0.99 

-  600. 

■  0.01 
-  2000. 

-  .015 

-  1000 

*  0. 
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The  value  of  v  a  1000  corresponds  to  a  filter 
cutoff  frequency  approximately  ten  times  the  largest 


observed  electron  bounce  frequency  in  the  ponderomotive 
potential  well  and  thus  effectively  models  an  infinite  band 
width  filter;  that  is,  there  is  no  frequency  discrimination. 

These  values  for  the  dimensionless  parameters  were 
chosen  to  simulate  an  FEL  with  the  following  physical 
parameters : 


Electron  Energy  50  MeV 

Electron  beam  current 

averaged  over  a  slippage  length  50  Amps 

Wiggler  Wavelength  1.25  cm 

Number  of  wiggler  wavelengths  400 

Signal  Wavelength  1  u 

Wiggler  Field  Amplitude  7  XG 

Nominal  Trapped  Electron  Energy  Loss  20  MeV 

Constant  Parameter  Length  7.5  cm 


A  constant  parameter  wiggler  of  fractional  length 
Vj  was  added  to  the  front  of  the  variable  wiggler  to 
enhance  the  linear  gain  per  pass.  Re  6  ■  .071  for  this 
case. 

In  Figure  6a,  the  normalized  integrated  pulse 

energy 

p  =  - ^  —■  j*  |  a  (v)  |  2  dv  (66) 
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is  plotted  as  a  function  of  the  pass  number  through  the 
wiggler.  It  may  be  noted  that  if  P  =  1,  the  energy 
lost  on  reflection  at  the  mirror  is  equal  to  the 
energy  extraction  from  the  electron  beam  assuming  all 
the  electrons  are  trapped  in  the  ponderomotive  potential 
well  formed  by  the  combined  field  of  the  pulse  and 
wiggler.  In  practice,  P  <  1  since  all  the  electrons  are 
not  trapped  and  additional  sinks  for  energy  loss  may  occur 
(for  example,  if  the  pulse  is  subjected  to  frequency 
discrimination;  see  the  next  section) . 

The  pulse  energy  P  increases  exponentially  at 
small  pulse  amplitudes.  After  nonlinear  saturation  of 
pulse  growth  occurs,  the  pulse  energy  P  oscillates  in  a 
limit  cycle  about  a  finite  mean. 

The  pulse  shape  at  saturation  is  displayed  in 
Figure  6b.  The  part  of  the  pulse  extending  from  v  ■  0 
to  v  *  1  interacts  with  the  electron  beam  in  the 
wiggler.  With  positive  phase  advance  8  >  0,  the  pulse 
is  shifted  to  the  left  on  each  successive  pass.  The 
front  part  of  the  pulse  (v  <  0)  shifted  out  of  the 
wiggler  no  longer  "sees”  the  electron  beam  and  is  damped 
by  reflection  energy  losses  at  the  mirrors .  Note  that  only 
the  front  of  the  pulse  is  amplified. 


The  Fourier  transform  of  the  pulse  amplitude  a(v) 


j”  eL  V  a  (v)  dv 


(67) 


is  plotted  as  a  function  of  <  in  Figure  6c.  For  a 
particular  value  of  k,  the  corresponding  frequency  w  of 
the  laboratory  frame  is 


to 

to 

s 


1  + 


k  L 
w  w 


(68) 


The  location  of  the  electrons  in  phase  space 


Y  -  r  dv, 


(69) 


is  plotted  in  Figure  6d  at  different  positions  v  in 
the  wiggler.  The  electron  energy  is  plotted  as  a  function 
of  the  phase  angle  modulo  2*.  If  the  electrons  are  freely 
accelerated, 


A 

Y 


|  rdv 


(70) 


where  y  is  the  initial  value  of  y,  and  is  set  to  zero 

o 

for  these  simulations.  If  the  electrons  are  trapped  at 
the  front  of  the  wiggler  and  remain  trapped  throughout 
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the  length  of  the  wiggler,  y  =  y Q,  and 


As  may  be  noted  from  Figure  6d,  the  electrons  are 
not  effectively  trapped  throughout  the  length  of  the 
wiggler  and  the  energy  transferred  to  the  pulse  at  satura¬ 
tion  is  small,  P  a  0.024. 

The  same  final  state  was  obtained  by  starting  the 
simulation  run  with  an  initial  square  wave  pulse  of 
finite  amplitude  large  enough  to  trap  the  begun  electrons 
(see  Figure  7).  The  back  of  the  pulse  erodes  away,  and 
sideband  frequencies  k~60  below  the  main  pulse  frequency 
(that  is,  zero  frequency  in  Figure  7)  grow  in  magnitude 
and  finally  dominate  the  spectrum.  The  sideband 
frequencies  are  of  the  order  of  the  electron  bounce 
frequency  ~|a|Js  ~  78  in  the  ponderomotive  potential  well. 
The  growth  of  sidebands  leads  to  a  loss  of  coherence, 
and  the  beam  electrons  are  not  effectively  trapped 
throughout  the  length  of  the  wiggler.  The  amplitude  of 
the  back  of  the  pulse  falls  to  a  low  level  compared  to 
the  peak  amplitude  at  the  front. 

This  suggests  that  any  tendency  of  the  pulse  to 
broaden  towards  the  back  during  growth  from  noise  will  be 
hindered  by  the  growth  of  sideband  frequencies.  Such 
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Figure  7a.  PsIm  Enargy  va.  Paee  Nunbar 


Figure  7b.  taplituda  i  va 


Figure  7c.  Complex  Amplitude  va.  Frequency  * 


Figure  7d.  Enargy  Extracted  va.  Phaaa  * 


figure  7g.  Energy  Extracted  va.  Phaae  * 


figure  7h.  Amplitude  va.  v 


sidebands  have  been  predicted  to  be  unstable  by  Kroll, 
et  al.2 

As  we  will  see  in  the  next  section,  the  suppression 
of  sideband  growth  by  passing  the  pulse  through  a  band 
pass  filter  results  in  a  considerable  improvement  in  the 
pulse  shape  and  the  total  pulse  energy  at  saturation. 

2 •  Frequency  Discrimination 

In  FELs  with  frequency  discrimination  where  the 
EM  pulse  is  passed  through  a  band  pass  filter  in  order  to 
attenuate  undesirable  sideband  frequencies,  the  band 
width  v  of  the  filter  appears  as  an  additional 
parameter.  Frequency  discrimination  by  a  band  pass 
filter  has  an  effect  proportional  to  l/v  which  tends 
to  retard  the  pulse  advance,  that  is,  it  tends  to  reduce 
the  effective  pulse  advance  produced  by  a  finite 
positive  value  of  B  .  To  avoid  the  EM  pulse  running 
either  ahead  of  or  lagging  behind  the  electron  beamlets, 
it  is  indicated  that  a  choice  of  v  *  1/8  is  desirable. 
This  was  confirmed  by  the  numerical  simulations.  A  small 
value  of  v  is  desirable  to  suppress  the  sidebands; 
however,  linear  gain  is  reduced  if  v  <  1/6  (see  Table  1). 
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TABLE  1 


u 


Variation  of  linear  gain  Re  5  with  v . 


r  -  0.99  n  -  600  T  ■  2000 


6  *  0.01 

Vi  * 

0.015 

Y  *  0  y 

th  “  ° 

v  160 

140 

120 

100  -  1/6 

90 

80 

Re  6  .017 

.0294 

.0367 

.0287 

.0143 

O 

. 

o 

1 

In  the  simulations  we  discuss  hereafter,  we  will 
consider  only  the  case  of  v  ■  1/6. 

With  v  *  1/6,  the  linear  gain  per  pass  for  a 
very  long  constant  parameter  wiggler  (  6/n  -*•  0)  is 

given  by^ 

Re  5 

z 

: .  i 


Re 


2(l-r+5) 


(nSz)!5 


=  2\ 


(72) 
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r 
*  * 


where  *th  -  0  and  Y  =*  0.  If  the  beam  and  wiggler 
parameters  are  kept  constant.  Re  [  B^/n*5]  increases 
from  zero  with  increasing  B,  reaches  a  maximum 

Re  [  BSj/n^j  =  .012  at  B/n*5  *  .042,  and  finally  decreases 

l' 

to  zero  at  B/n*  *  0.118.  The  theoretical  linear  gain 

curve  is  plotted  as  a  solid  line  in  Figure  8  while  the 

linear  gain  obtained  in  the  simulations  is  represented  by 

the  circled  points.  The  agreement  between  theory  and 

simulation  improves  as  B  is  made  smaller. 

If  the  wiggler  is  made  variable  (T  4  0),  the 

linear  gain  remains  close  to  that  of  the  constant  parameter 

wiggler  when  82r2/n  <<  1.  As  the  magnitude  of  the 

variation  increases.  Re  <S2  decreases  and  asymptotes  to 

zero.  The  theoretical  linear  gain  curve ^  for  a  very 

long  variable  wiggler  (B/h^  <<  1/  T  <<  1)  is  plotted  in 

Figure  9  as  a  solid  line  indicating  the  variation  of  6 

as  a  function  of  Br/n1*.  The  results  of  the  simulation 

are  represented  as  circled  points  for  values  of  B/n  2  ** 

-4 

5.8  x  10  ,  .01.  The  differences  between  the  observed 

linear  gain  and  the  theoretical  linear  gain  curve  for 
B/n**  *  .01  are  due  to  finite  length  effects  not  included 
in  the  theory. 

The  reduction  of  linear  gain  due  to  a  finite 
energy  spread  of  the  beam  becomes  significant  when 
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®  3  =  0.1  Simulations 


Figure 


8.  Linear  Gain-Constant  Parameter 

Wiggler  with  Frequency  Discrimination 


Figure  9.  Linear  Gain  -  Variable  Parameter  Wiggler 
with  Frequency  Discrimination 


T 


Yfch/r^  £  l.  This  may  be  seen  from  Figure  10  where 

Re  <S  is  plotted  as  a  function  of  7  ./T^  for  four 
2  th 

values  of  sr/n2.  The  energy  spread  of  the  beam  was 
modeled  by  a  Lorentzian  distribution  function. 

The  early  phase  of  the  evolution  of  the  EM  pulse 
growing  from  noise  levels  is  similar  to  that  discussed  in 
the  previous  section. 

For  a  variable  parameter  wiggler,  the  linear 
eigenmode  grows  in  amplitude,  the  peak  amplitude  moves  from 
the  back  to  the  front  of  the  pulse,  and  the  peak  amplitude 
at  the  front  continues  to  increase  to  levels  where 
electron  trapping  in  the  ponderomotive  potential  occurs 
(a  ~  T) .  Thereafter,  the  pulse  width  broadens  towards 
the  back  in  addition  to  increasing  in  amplitude.  The 
final  phase  of  the  evolution  depends  critically  on  the 
band  width  v  of  the  filter.  If  v  is  small,  the 
pulse  broadens  into  a  square  pulse  (approximately)  of 
width  v  -  1  and  the  electrons  trapped  in  the  pondero¬ 
motive  potential  well  remain  trapped  during  its  passage 
through  the  wiggler. 

The  results  of  a  simulation  run  for  the  same 
parameters  previously  shown  without  frequency  discrimina¬ 
tion,  showing  the  evolution  of  the  pulse  energy,  pulse 
shape,  and  pulse  spectrum  are  presented  in  Figure  11.  The 
FEL  parameters  are 
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Without 


Figure  11a.  Amplitude  va.  v 


Figure  Ilf.  Complex  Amplitude  va.  Frequency  < 


r  »  0.99 
n  -  600 

8  *  .007 

r  =  2000 

Vj  =  .015 
v  »  1/.007 

7th  "  °‘ 

In  contrast  to  the  simulation  results  of  the 
previous  section,  the  pulse  spectrum  remains  narrow  with 
a  large  peak  at  zero  frequency.  Sidebands  have  been 
suppressed  and  the  large  amplitude  saturated  pulse 
propagates  coherently.  Electron  trapping  is  very 
efficient.  From  Figure  11 j,  it  may  be  seen  that  —20% 
of  the  electrons  on  exit  from  the  wiggler  have  y  -  J  r dv  —  0 
and  therefore  these  electrons  have  been  freely  accelerated 
through  the  wiggler.  The  remaining  ~80%  of  the  electrons 
have  y  -  Jrdv  —  -2000,  which  implies  that  these  have  been 
trapped  throughout  the  length  of  the  wiggler. 

The  energy  extracted  is  about  80%  of  the  maximum 
possible  in  contrast  with  2%  for  the  same  parameters 
without  frequency  discrimination.  About  14%  of  the  energy 
extracted  from  the  beam  is  lost  in  the  frequency 
discriminators . 

Results  for  pulse  energy  and  the  fraction  of 
electrons  trapped  at  saturation  are  presented  in  Table  2. 
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TABLE  2 


r  »  .99 


v  -  1/6 


6 

r 

n 

vi 

P 

e4rn 

is 

ft 

(l-r) 7 

032 

1000. 

300. 

.022 

.49 

561. 

.7 

02 

1000. 

300. 

.022 

.57 

219. 

.75 

01 

1000. 

300. 

.022 

.68 

54. 

8 

.79 

005 

1000. 

300. 

.022 

.06 

13. 

7 

— 

022 

2000. 

600. 

.015 

.54 

530. 

.67 

007 

2000. 

600. 

.015 

.69 

53. 

7 

.80 

005 

2000. 

600. 

.015 

.06 

27. 

4 

_ 

The  total  pulse  energy  is  increased  by  increasing 

the  electron  current  h .  The  efficiency  of  energy  transfer 

is  improved  by  increasing  r. 

Further  increase  in  total  pulse  energy  may  be 

obtained  by  increasing  the  band  width  v  »  1/6  up  to 

some  critical  value  v  ... 

crit 

As  v  is  increased  above  the  critical  value 

v  ,  ,  the  pulse  shape  at  saturation  exhibits  a 
crit 
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progressive  depression  in  the  amplitude  at  the  back  of  the 
pulse.  Figure  12  displays  the  results  of  a  simulation 
run  with  FEL  parameters  similar  to  Figure  11  but  with 
v  increased  slightly  to  v  =  1  /6  s*  1/.005.  The  electrons 
are  no  longer  effectively  trapped  throughout  the  whole 
length  of  its  passage  through  the  wiggler,  and  the  energy 
transfer  from  beam  electrons  to  the  pulse  is  less 
efficient.  At  the  same  time,  the  spectrum  is  no  longer 
narrow.  Sideband  frequencies  of  the  order  of  the  electron 
bounce  frequency  begin  to  appear  and  eventually  dominate 
the  spectrum. 

The  deterioration  of  the  pulse  is  attributable  to 
the  growth  of  the  sidebands.  When  v  <  vcrit,  the  filter 
is  narrow  enough  to  suppress  the  growth  of  the  unstable 
sidebands  and  the  saturated  pulse  propagates  coherently 
as  a  finite  amplitude  square  pulse.  We  discuss  the  impli¬ 
cations  in  Section  4  below. 

Linear  gain  can  be  enhanced  significantly, 
particularly  for  large  B,  by  adding  a  constant  parameter 
section  to  the  front  of  the  variable  parameter  wiggler. 

The  constant  parameter  section  however  cannot  be  too  long 
since  electron  trapping  may  be  adversely  affected.  The 
linear  growth  phase  in  the  constant  parameter  section 
will  tend  to  saturate  at  (a)*5  vt  #  ir,  and  this  amplitude 
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(a)^  —  it/v  may  not  be  large  enough  to  begin  to  trap  the 
beam  electrons. 

In  Table  3,  we  present  the  results  of  a  series  of 
runs  in  which  the  effect  of  varying  the  fractional  length 
of  the  constant  parameter  wiggler  is  investigated.  It 

appears  that  the  length  of  the  constant  parameter  wiggler 

k 

should  be  limited  by  r  v  £  1  to  ensure  effective  electron 
trapping.  The  enhancement  of  gain  is  not  sufficient  to 
alter  our  previous  conclusions. 


TABLE  3 


r  =  .99  n  «  60.  r  -  200.  6  =*.05 


r1'2  v 

i 

6 

P 

o. 

.006 

.37 

.283 

.01 

.38 

.566 

.02 

.38 

in 

CO 

e 

.023 

.37 

1.13 

.031 

.02 

1.41 


034 


01 


3. 


Thermal  Effects 


The  presence  of  a  finite  thermal  energy  spread 

Y  k  in  the  electron  beam  affects  not  only  the  linear  gain, 

but  also  the  fraction  of  electrons  trapped  at  saturation. 

As  may  be  seen  from  Figure  10,  the  linear  gain  is 

_  % 

significantly  reduced  if  Y../r  >  1. 

"  tn 

In  order  to  trap  a  large  fraction  of  electrons, 

_  2 

Yfc^/2  should  be  less  than  the  well  depth  of  the  pondero- 

motive  potential  well.  The  effective  well  depth  at 

saturation,  estimated  from  energy  conservation,  is 

a  ~  [f.  nr/{l-r) 1^,  and  thus  the  limit  on  energy  spread  is 
o  t 


n 


% 


(l-r)T 


(73) 


1 

1 

I 


If  the  smallest  values  of  r\/(l-r)r,  as  determined 
by  Equation  (80)  in  the  next  section,  is  substituted  in 
Equation  (73),  the  limitation  on  energy  spread  may  be 
written 


(74) 
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Good  linear  gain  therefore  restricts  the  energy  spread  to 

about  half  that  which  would  be  deduced  from  looking  at  the 

requirement  for  effective  electron  trapping. 

Figure  13  summarizes  the  results  of  a  simulation 

using  an  electron  beam  with  Y  .  =  25.  The  other  FEL 

th 

parameters  are  identical  to  that  of  the  simulation 
presented  in  Figure  11. 

The  linear  gain  was  reduced  by  a  factor  of  3  from 

Re  6  *  0.029  =  0)  to  Re  6  as  0.01  ( Y.,  -  25.0).  Note 

tn  th 

that  in  these  runs  a  constant  parameter  section  of  length 

v  =  0.015  was  included  to  enhance  the  linear  gain. 

The  fraction  of  electrons  trapped  was  reduced  from 

f*.  »  0.8  (7*0)  to  f  ~  0.6  ( Y  *  25),  while  the 

^  th  t  th 

integrated  pulse  energy  was  reduced  by  75%  from 

P  -  0.69  (7*0)  to  P  -  0.52  ( Y  -  25). 

th  th 

These  results  support  the  conclusion  that  the 
limitation  on  beam  energy  spread  will  probably  be 
determined  more  by  the  need  for  linear  gain  than  by  that 
for  effective  electron  trapping. 

4.  Discussion 

We  may  interpret  the  results  of  Section  2  (7^*0) 

2 

in  terms  of  the  theory  of  sideband  instabilities,  which 
indicates  that  modes  of  frequency  u,  0.2  <  u>g  -  u>  < 
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Jj 

where  *  a  is  the  bounce  frequency,  should  be  unstable 
with  gains  several  times  the  signal  gain,  leading  to 
detrapping. 

To  suppress  sideband  frequencies  Aw  of  the  order 
of  a  quarter  of  the  electron  bounce  frequency  |  a  |  **  , 
the  reduction  in  the  pulse  energy  due  to  filtering,  i.e.. 


i;  |2 

r,  _  i.  1 

1  Aw  1 

i  ♦(£>*) 

[see  Equation  (32)]  should  exceed  the  relative  gain  due  to 
sideband  instabilities.  Thus 


>  2(1  -  r)  CG  -  1)  .  (75) 

where  G  is  the  ratio  of  gain  of  sideband  instability 
to  the  gain  of  the  pulse  amplitude.  Substituting  for 
the  amplitude  at  saturation  a2  -**  nrft/(l-r)  and  taking 
v  -  1/8: 


1* 


rn 


(1-  r) 


32  (G-  1) 

“FT” 


*76) 


Theoretical  estimates  indicate  3  «  2  to  3. 
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From  a  survey  o£  the  simulation  results  for  a  range 
of  FEL  parameters,  summarized  in  Table  2,  it  was  estimated 
that  for  stable  propagation. 


rn 


(1  -  r) 


S  50 


(77) 


where  the  numerical  coefficient  is  approximate. 

This  appears  in  good  agreement  with  the  above  theory. 

It  may  be  noted  that  large  3  (that  is  narrow 
band  width)  is  desireable  for  stable  propagation  in  the 
presence  of  the  sideband  instability.  However,  linear 
gain  becomes  small  if  B  is  too  large  (see  Figure  7). 
From  the  linear  eigenmode  analysis  of  a  very  long 
variable  parameter  wiggler,  the  linear  gain  Re  i2  is  a 
known  function  of  the  dimensionless  parameter  Br/n  , 

Re  6  =  fisr/n^).  Thus,  to  obtain  linear  gain,  it  is 
required  that 

^  <  (£)  *  ($)  • 


Equation  (77)  may  be  rewritten  as 

(i  -  r)r  ,  /  B  r  \4/3 


By  equating  the  right-hand  sides  of  Equations  (78)  and 
(79),  which  may  be  done  graphically,  we  determine  that  the 
inequalities  can  only  be  satisfied  ift 

ILUlL  <  o.o6  (80) 

n 

to  ensure  that  a  value  of  frequency  discrimination  may 
be  chosen  which  allows  both  for  linear  gain  and  stable 
propagation. 

From  the  definitions  of  r  and  n  (Equation  2^ 
and  Equation  28),  it  may  be  seen  that  the  above  inequality 
imposes  the  following  lower  limit  on  the  electron  beam 
current 


< 1  > 


1.42X105 

>  k  2  r  2 
w  o 


(1  -  r)  Ayr 


amps 


(81) 


where  is  the  beam  current  averaged  over  the  slippage 

distance  (uQ  ■  1),  the  wiggler  length  1^,  »  Tcsr02  is 
equal  to  the  "diffraction"  distance,  and  ft2/k2c2  »  1.  If 
the  beam  radius  is  restricted  to  £  *  so  that  the 

transverse  variation  of  the  pulse  amplitude  across  the  beam 
is  small,  the  minimum  current  criterion  is 
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6  r 

(i)  >  1.28  x  10  (1-r)  — f  amps 


(82) 


This  minimum  current  implies  a  minimum  circulating 

power  P  .  of 
min 


P  . 
mm 


(83) 


In  the  case  of  the  simulations  presented  in 
Figure  11, 


P  ,  »  84  GW 

min 

This  value  of  P  is  about  two  times  the 
min 

estimate  of  minimum  circulating  power  calculated  by  Kroll, 

2 

et  al.  for  Ay  r/yr  ■  0.4.  Their  calculation  takes  account 
of  the  constraints  imposed  on  pulse  amplitude  variations 
transverse  to  (kwr0  ■  1/3)  and  along  (Lw  *  ksrQ2) 
direction  of  propagation,  but  does  not  include  the  need  for 
good  linear  gain  and  suppression  of  unstable  sideband 
frequencies. 
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5.  Summary 

We  have  developed  a  1-D  code  for  the  simulation  of 
FELs  following  the  self-consistent  interaction  of  electrons 
and  the  EM  wave  in  a  wiggler,  while  reflecting  and 
optically  filtering  the  EM  pulse  for  multiple  passes 
through  the  system.  In  this  way,  both  pulse  formation 
and  the  steady  state  may  be  studied.  We  have  included 
provisions  for  both  constant  and  variable  wiggler s  and 
for  thermal  energy  spread  in  the  entering  electron 
pulse.  In  this  paper  we  concentrate  on  the  results 
obtained  with  electron  pulses  of  zero  length 
(6 -function)  which  reduces  the  simulation  time  required. 
This  should  be  a  good  representation  for  pulses  such  that 

u  <1. 

0 

Our  principal  results  are: 

(1)  We  have  verified  analytic  results  previously 
obtained  for  linear  pass  to  pass  growth  of  the  pulse  at 
low  amplitude. 

(2)  We  find  that  for  high  saturation  systems 
optical  discrimination  is  required  in  order  to  prevent 
nonlinear  signal  breakup.  Qualitative  agreement  with  the 
theory  of  sideband  instabilities  is  obtained  for  the  width 
of  the  required  frequency  filter. 

(3)  The  further  requirement  that  the  linear  phase 


of  growth  saturate  at  a  high  enough  level  to  ensure 
trapping  seems  to  be  satisfied  except  for  systems  which 


seek  to  enhance  linear  gain  by  use  of  a  long  constant 
parameter  section. 

(4)  With  proper  choice  of  frequency  discrimination 
(and  for  an  ultrashort  electron  pulse)  we  have  demonstrated 
very  high  extraction  systems  (30%  efficiency  at  lv) 

which  grows  from  noise  to  a  stable  steady  state.  We 
derive  and  numerically  verify  a  criterion 

g  (1  -  r)  Ayr 

<I>  >  1.28  x  10  - 5 -  Amps 

Yr 

for  the  minimum  current  required  for  such  a  system. 

(5)  The  limitation  on  beam  thermal  energy  spread 

Yfch  is  determined  more  by  the  need  for  good  linear  gain 

(Yth/rJs<  1)  than  by  that  for  effective  electron  trapping, 

(7/  rh  <  2). 

cn 
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